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The Brownian motion of a quantum particle in a harmonic confining potential and coupled to a 
harmonic quantum thermal bath is exactly solvable. Though this system presents at large temper- 
ature a pedagogic example to explain the laws of thermodynamics, it is shown that at low enough 
temperatures the stationary state is non-Gibbsian due to an entanglement with the bath. In physical 
terms, this happens when the cloud of bath modes around the particle starts to play a non-trivial 
role, namely when the bath temperature T is smaller than the coupling energy. Indeed, equilib- 
rium thermodynamics of the total system, particle plus bath, does not imply standard equilibrium 
thermodynamics for the particle itself at low T. 

Various formulations of the second law are found to be invalid at low T. First, the Clausius 
inequality can be violated, because heat can be extracted from the zero point energy of the cloud of 
bath modes. Second, when the width of the confining potential is suddenly changed, there occurs 
a relaxation to equilibrium during which the entropy production is partly negative. In this process 
the energy put on the particle does not relax monotonously, but oscillates between particle and 
bath, even in the limit of strong damping. Third, for non-adiabatic changes of system parameters 
the rate of energy dissipation can be negative, and, out of equilibrium, cyclic processes are possible 
which extract work from the bath. Conditions are put forward under which perpetuum mobile of 
the second kind, having one or several work extraction cycles, enter the realm of condensed matter 
physics. Fourth, it follows that the equivalence between different formulations of the second law 
(e.g. those by Clausius and Thomson) can be violated at low temperatures. 

These effects are the consequence of quantum entanglement in the presence of the slightly off- 
equilibrium nature of the thermal bath, and become important when the characteristic quantum time 
scale h/ksT is larger than or comparable to other timescale of the system. They show that there 
is no general consensus between standard thermodynamics and quantum mechanics. The known 
agreements occur only due to the weak coupling limit, which does not pertain to low temperatures. 

Experimental setups for testing the effects are discussed. 

PACS: 05.70Ln, OS.lOGg, 05.40-a 



1. INTRODUCTION 



The faith in the laws of thermodynamics has been strengthened time and again because numerous counterarguments 
and perpetuum mobile setups failed. It was summarized in the classical statement of Arthur Eddington in 1948 [Q: 
"The law that entropy always increases - the second law of thermodynamics - holds, I think, the supreme position 
among the laws of Nature. If someone points out to you that your pet theory of the universe is in disagreement with 
Maxwell's equations - then so much the worse for Maxwell's equations. If it is found to be contradicted by observation, 
well, these experimentalists do bungle things sometimes. But if your theory is found to be against the second law of 
thermodynamics I can give you no hope; there is nothing for it but collaps in deepest humiliation. Nevertheless, what 
we intend to do in this paper is to show that several formulations of the second law need not apply to systems coupled 
to a bath in the quantum regime. This paves the way for a new, modest definition of the most despised objects of 
modern physics, perpetuum mobile of the second kind. We shall propose realizations wherein they can make a few 
or even many cycles, though not infinitely many. A short version of the material appeared already ||^, which was 
discussed in the scientific press 

The laws of equilibrium thermodynamics apply both to (quasi) closed quantum and classical systems, and to open 
classical subsystems |^ . This can all be traced back to the general character of the Gibbs distribution that describes 
the equilibrium state. The same laws are believed to apply as well to open quantum subsystems. Our aim will be to 
show that, though this belief is proper for weak coupling, it is not justified for non-weak coupling between system 
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and bath. Non-weak coupling means physically that a cloud of bath modes has been formed around the particle, 
which we shall still consider as part of the bath. Although one could defend the opinion that — as an example — for 
Kondo-type problems there is no need or no sense to describe the impurity spin ("the subsystem") separately from 
the Kondo-cloud of conduction electrons ( "the bath" ) , such a viewpoint is not obvious for an array of mesoscopic 
Josephson junctions, where the "bath" is the electro-magnetic field, to which they may be strongly coupled. Since 
that regime shows up in many systems at low enough temperature even for a small but fixed coupling, there is an 
important case to make. 

When deriving the Gibbs distribution for a (quasi)closed system ||,^, the underlying statistical assumptions do 
not depend much on the quantum or classical nature. An open subsystem is usually considered as being in contact 
with an initially equilibrium bath. Under some general statistical conditions concerning the bath ||-|,|n|,|l|l , which 
are again the same for the quantum and classical cases, one derives a Langevin equation. The general behavior of 
the classical Langevin equation is well known |^ : The stationary distribution is Gibbsian, and, at least for the white 
noise case, the equivalent Fokker-Planck equation is the main tool to describe the non-equilibrium statistical theory 

£!■ 

Much less is known about the quantum Langevin equation [^,0jp-|lJJl^ , first proposed by Senitzky in the weakly- 
damped (weakly coupled) case, and in a more general form by Ford, Kac and Mazur Its stationary distribution 
has been obtained for the harmonic potential, in which case it depends explicitly on the coupling constant, becoming 
Gibbsian only in the limit of weak coupling |l4|Jl^ , p^ . 

Let us recall that the situation with a particle interacting with an equilibrium bath is known as Brownian motion, 
and the particle as well is referred to as a Brownian particle. As one of the paradigms of non-equilibrium statistical 
mechanics, the theory of Brownian motion has numerous applications in condensed matter physics |]6|, pT| , p"3|JT^Jl^ , p^ , 
atomic physics [p|Jl9|], quantum optics and chemistry [^,^. It is believed that some of those practical realizations 
can be understood as being weakly-coupled with their thermal baths, and then standard methods of statistical physics 



can be applied [£1 12|. However, there are nowadays well-known experimental situations, which are essentially far 



from the weak-coupling regime. The main example is the case of weak links between superconductive regions, so- 
called Josephson junctions, in their overdamped regime [ p^JT?! , where the relevant ranges of parameters were achieved 
already twenty years ago. Even in quantum optics, which has often been satisfactorily described by weak-coupling 
theories [ p^ , there are recent experiments showing the necessity of moderate and strong coupling approaches (see, 
e.g., 1^^). In NMR experiments, on the other hand, very weak coupling occurs, but it leads to a T2 relaxation time 
of transverse (non-classical) correlations, that can reach up to minutes. During that time related non-thermodynamic 
effects can occur. 

The cause of the crucial differences between quantum and classical Brownian motion lies in quantum entanglement: 
A complete description in terms of a wave function is possible only for a closed system; subsystems are necessarily in 
a mixed state. Since the quantum Gibbs distribution tends to the pure vacuum state for T ^ 0, it is in that limit 
not an adequate candidate for the description of the quantum subsystem non-weakly interacting with its thermal 
bath. Where typically researchers have just guessed that thermodynamics would apply anyhow, we shall provide the 
opposite answer by analyzing exactly solvable models. 

Another important property is embodied in the structure of the quantum Langevin equation : As predicted by 
the quantum fluctuation-dissipation theorem, but in contrast to the classical case, the time-scales of fluctuations and 
dissipation are different, and even in the limit of instantaneous friction, the noise does not become white, but has 
the characteristic quantum timescale Ti/ksT . During this period the noise has a memory and thus has the possibility 
to cause non-statistical but "mechanical" or "quenched" , in any case non-thermodynamic, behavior of the system on 
which it acts 1^,0]. The physical meaning of the classical fluctuation-dissipation theorem thus is stronger, since it 
insists on equal timescales of friction and noise, which induces the standard thermodynamic behavior. 

The new properties become non-effective if the interaction with the thermal bath is weak or if temperature is high, 
in which case the subsystem relaxes to its own quantum or classical Gibbs distribution. Both these cases will be 
referred to as the Gibbsian limit. We recall that its existence is typically not connected with very low temperatures, 
since even for a small but generic coupling between the system and the thermal bath, one always goes out of the weak 
coupling limit by taking the temperature low enough. 

It is a crucial fact that a non-Gibbsian stationary state implies the inadequacy of equilibrium thermodynamics. 
In the present paper we propose a generalized thermodynamical description of a quantum Brownian particle in 
a harmonic confining potential. This description is self-consistent, and docs not use a priori the concepts of the 
equilibrium (Gibbsian) statistical thermodynamics. Instead we employ the methods of nonequilibrium statistical 
thermodynamics developed recently for glasses and applied to black holes, where a separation of timescales allows for 
a two-temperature approach [p3|-p6[ . 

The universal character of equilibrium thermodynamics led to the general expectation that in one way or another. 
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thermodynamics will be applicable to the full quantum domain. A somewhat stronger point was expressed by Landau 
and Lifshitz jsj, namely that the proper formulation of equilibrium thermodynamics must be based on quantum 
mechanics. For the strongly coupling quantum situation one might, however, not be convinced. Let us give three 
principle arguments that question standard wisdom: 

1) The bath modes are correlated during the quantum timescale h/T, even when damping is instantaneous. When 
this timescale is larger than the (largest) relaxation time of the system, the bath acts more like a quenched disorder 
than as a white noise. Thus the standard condition for going from a Langevin equation to a Gibbs distribution is not 
fulfilled and new behavior should be expected. 

2) Assume that the overall system (the Brownian particle plus the bath) is in equilibrium at a low temperature. One 
of the many formulations of the second law states that no heat can be extracted from the bath. This just follows 
from the Clausius inequality: dQ < TAS (here dQ is the heat flowing from the bath to the subsystem, the Brownian 
particle, while S is the entropy of the subsystem) for T — > 0. A naive argument in support of this statement will 
be that the bath is close to its ground state, and therefore its energy cannot be lowered. However, this argument is 
invalidated by quantum entanglement: Since the bath is not an isolated system itself, it cannot be in a definite energy 
eigenstate. In particular, it cannot be in the ground state, and its energy fluctuates even for T 0. 

3) If a closed quantum system is its ground state, the only change can be to do work on it. Now suppose that this 
system consists of a subsystem coupled to a bath, and that the external coupling connects to the subsystem, and not 
to the bath. Then typically the action of doing work will reshuffie also the separate energies of the subsystem and the 
bath. As the direction of the exchange depends on the question whether externally work was added or subtracted, in 
one of the cases the subsystem receives energy from the bath. Since this comes from the unobserved bath modes, it 
cannot be identified as work, and must be identified as heat extracted from the bath, even at T = 0. 

Because these arguments question standard wisdom, the only way to investigate the practical situation is to start 
from first principles, namely from standard quantum mechanics. This is the general strategy of the present paper. 

Statements on violations of certain formulations of the second law in the quantum micro world already appeared 
in literature. Capek and his coworkers noticed such effects in certain biophysically inspired models, and Nikulov 
pof reported on violations of the second law in mesoscopic superconductivity. The latter author bases his view on 
results for permanent currents in inhomogeneous superconducting rings |Q . 

Since the subject of violating the second law and introducing perpetuum mobile has such a notorious history, 
new works in this field should be as convincing as possible. Therefore we present now a somewhat extensive, but 
self-contained exact analysis that leads to our conclusions, partly already presented in Our methods are exact 
since the case of a quantum particle in a harmonic confining potential and coupled to a bath of harmonic oscillators 
with bilinear couplings is exactly solvable. Notice that in a previous work Q we also entered the discussion of the 
approximate solution for a weakly anharmonic force; this will not be touched here. 

Our paper is organized as follows. In section 2 we recall the derivation of the quantum Langevin equation. In 
section 3 we solve the statics of the total system (subsystem plus bath) by elementary techniques. In section 4 we 
show that the thermodynamics of adiabatic changes can be described through two effective temperatures, and that 
analogies with the usual two temperature thermodynamics can be stated much further: The generalized relations 
will have the effective form of the first and second law. Next we first present details on the violation of the Clausius 
inequality. In section 5 we consider the dynamics of the system for the case where the initial state is Gibbsian, 
and for which the spring constant is slightly modified at time zero. In section 6 we use those results to derive the 
energy relaxation and the entropy production. In section 7 we consider work done on the system for that situation 
of an instantaneously changed spring constant, and we also consider work for smooth, cyclic changes. In section 8 
we mention a number of experiments where our results have applications. In section 9 we view those results from 
the foundations of thermodynamics and the definition of perpetuum mobile. Finally in section 10 we close with a 
discussion. 



2. THE QUANTUM LANGEVIN EQUATION 

It is standard wisdom that the analysis of Brownian motion of non-interacting particles may be restricted to a 
single Brownian particle. This insight goes back to Szilard in his somewhat related work on Maxwell's demon 
for a translation see | |3^ . In our analysis we shall also make this simplification, but insist that the bath has many 
degrees of freedom. Therefore it is equally equipped to couple to a gas of N non-interacting Brownian particles, and 
our results for energy, work, heat, entropy, etc. must just be multiplied by N . Because of this, our results yield, 
without say, the intrinsic variables of a large Brownian system in its thermodynamic limit. 



3 



1. The Hamiltonian 



The quantum Langevin equation is derived from the exact Hamihonian description of a subsystem (Brownian 
particle) and a thermal bath, by tracing out the degrees of freedom of the bath. The influence of the particle on the 
bath is assumed to be sufficiently small. Thus, only the linear modes of the bath are excited, and the interaction of 
the particle with the bath is assumed to be linear. To be as pedagogic as possible, we first take a definite model for 
the bath, namely a collection of harmonic oscillators; later we will relax this assumption. For the total Hamiltonian 
we thus assume Wu 



'Htot = Ti. + TLb + 'Hi 
H = IC{p) + V(x), 



2m' 



(2.1) 
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where H is the Hamiltonian of the particle, consisting of the kinetic and potential energies /C and V, Hb is the 
Hamiltonian of the bath, and Tij is the interaction Hamiltonian. p, pi, x, Xi are the momenta and coordinate 
operators of the particle and the linear modes of the bath. V{x) is the confining potential of the particle, and m and 
rrii are the corresponding masses. 

Notice that our total system is closed and its energy is conserved, except for the periods where work is done on it 
by externally changing a system parameter such as m or a. When we later on take as initial density matrix the Gibbs 
distribution exjp{— (3Ti.tot) / Z , this still refers to our closed system; in particular it is not part of a larger thermodynamic 
system, with which heat exchange would be possible . 

The Si-terms of Tis + 7i/ form a complete square, since Tij includes a self-interaction term proportional to . This 
guarantees that the total Hamiltonian Titot will be positive definite. In certain physical situations (e.g. electromagnetic 
interaction), such a term is generated automatically by the coupling Indeed, under a canonical transformation: 



P^ 



-XiJUiUJi 



(2.2) 



Htot becomes 



'^tnt = 



2m 



E 



2m 



■{Pt 



2 , mjUJ^ 2 



(2.3) 



which corresponds to the minimal coupling (or subtraction) scheme. 

In other situations (such as certain cases in nuclear and atomic physics, see pO| for more details) the self-interaction 
term is absent, and one has H 



p 

2m 



V{x) 



•Hi 



CiXi 



(2.4) 



while Hb is unchanged. For a harmonic system one will have V{x) — ■^bx'^. In general, the potential energy will have 
a minimum only if 6 = V" (0) is large enough. This happens when the combination 



■niiUJf 



(2.5) 



is positive. In the case we shall consider below, with the Ci given by ( 2.19| ), this sum can be evaluated, leading to 
a = b — 7r, where 7 is the coupling strength (damping constant) and F is a large cut-off frequency. So this system can 
be mapped on previous one provided we define V{x) — V{x) + ^7ra;^. In doing so we identify with Hamiltonian of 
the subsystem the expression 7Y = /C + V of (2.1), rather than /C + V, and with Hi the combination Hi + ^jTx'^. To 
give some motivation for this, let us notice that we shall take F large and 7 finite. In that case the expectation value 
(V) is large, proportional to 7F, but this is almost fully compensated by an opposite term arising from (Hi). These 
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cancellations have been accounted for by the mapping to 7i and 7i/, leaving at most a InF divergence for large F, 
which actually arises at small temperatures from {K) . We shall come back to pro's and contra's of this identification 
in section 4. 2e, where we notice that it is already needed to obtain the standard thermodynamics at very large T. 



Some word of nomenclature is called for. The case of a harmonic potential V{x) = ^ax^ is often called "the linear 
case" in literature, of course referring to its linear force. The expressions "linear potential" and "non-linear potential" , 
that are sometimes found in literature, are misnomina, and will be avoided by us. 



2. Derivation of general quantum Langevin equation 



The operator equations of motion for the bath modes read 



. _ 1 

Xi — Pi 

rui 



Pi = -XiHTliUJi + CiX 

After introducing the creation and annihilation operators by 



V zniiOJi V 2 



we can write Eqs. (2.6, 2.7) in the form 



(2.6) 
(2.7) 

(2.8) 



i{t) = -iujiai{t) + i\ 



2hmiUj. 



■x{t). 



(2.9) 



This equation is solved readily: 



a.(i) =e-*"'*a,(0)+z^ 



2hm.iUJi Jq 



dse-*"'(*-"'a;(s), 



(2.10) 



yielding 



Xi (t) = Xi (0) cos Wit H sm Uit H 

niiOJi niiUJi Jq 



ds sin.uJi{t — s)x{s) 



Pi{t) = —miUJiXi{0) sinujit + pi{0) cos ujit + Ci / ds cos uji{t — s)x{s) 



(2.11) 
(2.12) 



The Heisenberg equations of motion for the Brownian particle read 

1 

x = —p, 
m 

p — ~ax + CiXi — x 



(2.13) 
(2.14) 



Combined with Eq. (2.11) the last equation becomes 



mx{t) = -ax{t) + rjit) - [ dt'j{t - t')x{t') - ^Yx(t) 

Jo 



where 



(2.15) 



ilit) = Vc,[a;,(0)cosc^,t + ^sinw.i] = V W 7r^[al(0)e^"'* + a,(0)e-*'^'*]. 



■niiUJi 



2m,w,: 



(2.16) 



(2.17) 
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are the noise related to the unperturbed bath, and the friction kernel, respectively. A partial integration brings 



mx{t) ^ -ax{t) - x{0)j{t) - [ ds-f{t - s)x{s) + r]{t), (2.18) 

Jo 

Notice that in this derivation the back-reaction of the bath on the particle has been taken into account in an exact 



manner. It is described by the integrals in eqs. ( 2.11 ), ( 2.12 ), and brings the damping terms x(t)j(0) — x{0)^{t) 
/(J ds"f{t — s)x{s) in going from eq. ( ^.14 ) to ( ^.18 ). 



3. Drude-Ullersma spectrum 

For some, but not all, of our applications it is benefitable to consider a fully explicit case for the bath. The bath is 
assumed to have uniformly spaced modes 

uj^=iA 1 = 1,2,3, ••• (2.19) 
and for the couplings we choose the Drudc-UUersma spectrum |TTp|] 



Here F is the characteristic Debye cutoff frequency of the bath, and 7 stands for the coupling constant; it has dimension 
kg/s. Our parameter 7, related to another one (7' = j/m) sometimes employed, see, e.g., ]lT[ , allows to consider 
changes in the effective mass m at fixed coupling to the bath. 

The thermodynamic limit for the bath is taken by sending A — > 0, which induces relaxational behavior. As usual, 
the "Heisenberg" timescale 1/A will be extremely large, implying that in the remaining approach the limit of "large 
times" always means the quasi-stationary non-equilibrium state where time is still much less than 1/A. In the limit 
A — each coupling Ci ~ -s/A is very weak. The fact that the bath has many modes nevertheless induces its non- 
trivial influence. At finite but small A the system would have an initial relaxational behavior, which at times of order 
1/A is changed in a recurrent behavior. 

It is customary to define the spectral density 

i 

It has the Ohmic behavior J "fuj for cj <C F, and 7 is called the interaction strength or damping constant. As J(w) 
is cut off at the "Debye" frequency F, it is called a quasi-Ohmic spectrum. 

For many applications only the spectral density needs to be defined. However, for our further calculations it is 



advantageous to stick to the fully spec ified linear bath, with its frequencies ( 2.19 ) and couplings ( 2.20 ). It can then 



be shown that the friction kernel (2.17) becomes 



00 



7(i) = — / duj / coscjt = 7Fe~^l*l (2.22) 



-r i 

It is non-local in time, but on timescales much larger than 1/F it may be replaced by 7(5+ (t). 

Finally we wish to mention that there are alternative ways to derive the quantum Langevin equation pl[ |, since 
many of its properties are rigidly determined by general statements like the quantum fluctuation-dissipation theorem 
Nevertheless, we choose to focus on concrete models, because they show in detail how the quantum Langevin 
equation arises from first principles, and thus are better suited for pedagogical purposes. 

3. GIBBSIAN STATE FOR A HARMONIC CONFINING POTENTIAL 

The case of an oscillator subject to a parabolic confining potential is a celebrated exactly solvable problem in 
quantum mechanics, 
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n{p,x)^JC{p)+V{x) = ^ + l-ax^ (3.1) 

2m 2 

The eigenfrequency is already know from the classical treatment, 

(3.2) 

When this oscillator is coupled to an oscillator bath with bilinear coupling, as was done in Eq. (^.l|), the problem 
remains exactly solvable. It is in the true sense "the harmonic oscillator model" for quantum Brownian motion. 

It is well known that, besides its direct physical meaning, the harmonic oscillator can be interpreted as an LC 
circuit ||l^. Then x may correspond to the charge Q on a capacitor, 1/a to its capacitance C, m to an inductance 
L, p to a flux LQ, 7 to a resistance R, and r]{t) to a random electro-motoric force. Although we will not use this 
language explicitly, it is useful to keep it in mind, especially when considering variations of parameters. Indeed, in 
this setup there should be nothing very difflcult in varying L or C, in our notation m and a. 

The theory of the dissipative harmonic oscillator is considered in many works (see [|ll| and Rcfs. there, as well 
as a recent work for the driven case p9|]). We will now be primarily interested in thcrmodynamical aspects of this 
problem. 




1. Shift of the bath frequencies due to coupling with the central particle 



In Fourier space the equation of motion of the particle may be written as 

^2 



(a — muj'^ + — = CiXi, 



and for the bath 

{—niiUj'^ + ■miLjf)xi — CiX 
From these relations one derives a condition for the cigcnfrcquencies Vk 



(3.3) 



(3.4) 



^22 \ ^ 

V = u } 

m — ^ 



^ mm,w2(w2 _ j^2) 



E 



A 



(c.2+r2)(c.2-^2) 



A 



El 



A 



(3.5) 
(3.6) 



where we inserted the definition ( 2.20| ) of the c^. The first sum may be replaced by an integral, while the second can 
be carried out exactly, 



El 



A 



A 



= lim 

UJi ^ V LOi + V' N^ca 



■^(7V-^)-^(l-^)-V(iV+^)+^(l + ^)] =^-^cot^ (3.7) 



where ipi^) — dlnr(z)/dz is the di-Gamma function and we used 

4>{z + 1) = V'(^) H 7 "0(1 ^ z) ~ '0(^) + TTCOtTTZ 

Z 

The eigenfrequencies i^k of the coupled system thus follow as the roots of 



TTV A 

cot — — 

A TTl/ 



(3.8) 



(3.9) 



The transcendental equation has no solution for Q < u < l^. For 1/ > A there is one solution in each period of the 
cotangens , except for the period that contains the point v — ujq = ^ a/m, where there occur either three solutions or 
one. One can then check that in the limit of vanishing coupling 7—^0, there occur the modes LVi = lA (z = 1, 2, • • •), 
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and uiQ. Notice, however, that this behavior only pertains in the regime of infinitesimal coupling 7 < 7c with jc ^ mA. 
For 7 > 7c, however, the interval containing uo has only one solution, so ujq is lost as a separate mode, its influence 
being taken by a shift of neighboring modes. 

For finite 7 the solution of (|3.9|) shows that the bath modes LVk '3> A now get shifted to 



where 



z/fc = fcA - -(t){kA)A = ujk - -(j}{ujk)A (3.10) 



(iy)=arctan- „ — --ttt — ^ (3-11) 



Here the definition of the arctan is such that <^ goes monotonously from (/)(0) = to 0(oo) = tt. We shall need 

7r2z. 



sin0(z^) 



7J/ 



{(a — rav^Y' + 72i/2}i/2 
where the approximation holds for large F. 



(3.12) 



2. The Gibbsian state of the particle and its bath 

The steps of last section allow to calculate the Gibbs free energy of the total system, 

/?i^tot(T,7) =^ln2sinhi/3?ii/fe (3.13) 
fc 

For small A one may use the identity 

E^w^^/ ^.M.,^^-l ^.^AM^l 



00 



MU-'-^]Ai.)+OiA) (3.14) 



k=l 

and one gets 

/3Ftot(T,7) = /?FB(r, 7 = 0) + /3i^p(a, 7, F,m,r) (3.15) 

where the first term is the free energy of the bath in absence of the particle. Neglecting its divergent zero point energy 
one gets 

PFb{T,j^ °) = a i, '^'"^''^^ " '""''"^ = ^^-^^^ 
It is of order 1/A, showing the extensivity of the bath, and implies the energy 



Ub{T,j = 0) = —T' (3.17) 



and the linear specific heat and entropy 



Cb{T,^ = 0) = Sb{T,j^O) = ^ (3.18) 
The free energy shift due to the central particle, its coupling to the bath and the resulting disturbance of the bath. 



follows from Eq. (3.14) as 

PF.p = - [ di^ ln[2sinh — (3-19) 

TT Jq 2 dv 

r A 1 [9 • 1. 1 flfi 1 ar2 + (mr2+7F-a)^2^3^^4 

= / di/ In 2 smh -/3H 77 2\( 2 , r2^ , r 212 , 2r4 2 ^-^O) 

n Jq 2 [(a — mj/^)(i/'' + 1 ^) + 7I j^^J^ + 7^1 

A useful identity is 

11/"°° 1 1 

/3Fp = In 2 smh{-(3hvi) + -/3h / [0{v - vi) - -(/)(i')]cotanh (-Phv) (3.21) 

2 2 Jq tt 2 

where ui is arbitrary and 9 is the Heaviside step function. 
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a. Intermezzo: the characteristic frequencies of the damped oscillator 



The present model for a damped harmonic osciUator has three characteristic frequencies, that do not depend on 
temperature. They just foUow from the Uncar equations of motion, and thus have the same value at high and low 
temperatures. 



The denominator in (3.20) is a fourth order polynomial in . It decomposes as 

to2(j,2 ^ r2)P3(ity)P3(-ij/) (3.22) 

where 



„3 -n„2 t .2 I , ,2-n _ /„ t^\i„2 , , 2\ , 7^ 



P3(s) =s'' - Ts' + {uof. + ^)s - ujiV ={s- T){s' + ui.) + ^ s (3.23) 

TO TO 

The roots s = ti'1,2,3 of ^3(5) satisfy the relations 

U1+LJ2+ = r (3.24) 
UJ1UJ2 + Lo-i^^s + ^3^1 = ^ (3.25) 

TO 

ujiuj2i^3, — Ttiig (3.26) 

Two different situations can arise: Either all three roots are real (this is the case in the overdamped regime), or, 
in the underdamped regime, two of them are complex conjugate: = W2, whereas 0^3 is real. In both cases one has 
Re Ci^i,2,3 > 0, which indicates that with time the particle relaxes toward a stationary state. 

For small 7 these roots read: 

+ ■ ,7 r , ^7^2 r2(r±»^o) 

'^^-^"TO(r2+c.2)-l;;^j (r2+a;2)3 ^3.28) 



On the other hand, for a large F one gets 



7 / / 4aTO \ 1 / 7 \ 



1 — 2amH'^ 
1± ^ 



(3.29) 



^3 = r - ^ - - (^) (3.30) 

TO I V TO / 

We shall only need them to leading order in 1/F 

7(1 -u;) 2a 7(1 + w) 7 

^1 = — 5 = .-I . V = — ^ , t^3 = r (3.31) 

2to 7(1 + w) 2m m 



where we denoted 



Later on we shall need 



am , 

e=^, w = Vl - 4e (3.32) 

7 



dw dw 1 — w dwi,2 7(1 — w ) da;i^2 , 7 



a— — = to- — = , a- 



da dTO 2w ' da 4toi(; ' dTO 4TO?i; ^ ~'~ ^ ^ 

For overdamping (e < i) w is real positive. Our interest is in particular the strong damping regime 7^ 3> am, 
where 

^1 = -, c^2 = — (1 5-), t^3 = r (3.34) 

7 TO 7^ TO 
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and the approximations hold to leading order in e. 

Already in the classical regime our system has three characteristic relaxation times: for the coordinate, for the 
momentum and for the noise. For large T and 7 they are well separated 



1 



7 



LUi a 

In the quantum regime the quantum time-scale 



1 m 

UJ2 7 



1 
f 



(3.35) 



T 



is comparable to or larger than t^, inducing quantum coherence effects of the noise and thus new physics. 
In case of underdamping, (e > i) one has w = iw, with 



w 



^1,2 



7(1 =F iw) 
2m 



This leads to the renormalized oscillation time tq and the damping time 



1 



To = 



2m 

1 



(3.36) 



(3.37) 



(3.38) 



Since differs from Tp by a factor of order unity, we may skip the latter and use taup and tq as the relevant 
timescales in the underdamped regime. 

It is worth to mention th at the weak-coupling limit commutes with the quasi-Ohmic limit, in the sense that taking 
large F i n Eq s. ( 3.27 , 3.28 ) we get the same main term and at least the first correction as having taken small 7 limit 
in Eqs. (||3|, |3.30|). 



h. Continuing the main argument for the Gihhsian state 



In order to calculate the free energy (3.20), we shall first determine the following integral 

dv V coth(iai^) 



We can write 



/(a,AS)= / 
Jo 



/(a, AS) 



Aw 



(A2 + I.2)(52+j,2)^ 



Aw 



(A2 +zy2)(52 + j,2) (e''" - l)(yl2 +ty2)(52 _^ j,2) 



\w 



\w 



B^-A-^ 

where we used the known formula 



(e'^''-l)(A2 + z.2) (e'^'' - 1)(S2 -h j/2) 

TT 1 

' aAB{A + B)' 



tdt 



In z — ipiz) 



(exp(27rt) - l)(i2 + z2) ' 2z 



By integration we obtain 

J{a,A,B) = ^ ln(2sinh-ai.) i ^2^^2 



B^ + i 

1, ^faA\ I aA 1, ^faB\ 1 , aB 

-— InF -In — - H InF H In — - 

A \2n 2A 47r2 B \2tz 2B Att^ 



(3.39) 



(3.40) 



(3.41) 



(3.42) 



In terms of the roots uji, we may write Eq. (3.1C) as 
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(y) = arctan \- arctan ■ 



UJ2 



V V 

arctan arctan — 

LJ3 r 



(3.43) 



The derivation follows immediate ly aft er using ( 3.22 ) with ^3(5) = (s — wi)(s — W2)(s — W3) and expressing the arctan in 
logarithm's. The integral in Eq. ( 3.20 ) can no w be done by adding to i>{v) a ter m \{T — ld\ — u}i^lo^)Ild/^ arctan(z^/w4), 
which vanishes for any CJ4 on account of Eq. ( [3.24 ), and then using Eq. ( 3.42 ) with a = hf3. 
This finally brings the shift of the free energy due the presence of the particle 



InE 



InE 



V 2tt 



InE 



V 2tt 



V 27r y (27r)2 



(3.44) 



This is just equal to — hiZ' with Z' calculated in Eq. (4.20) of Grabert et al. These authors did not point at 

the physical role of their Z' . Here we see it is the part of the partition sum of the total system related to the central 
particle and its coupling to the bath with its linear unperturbed spectrum cj^. = fc A. We nevertheless expect that the 
statics and the dynamics hold for more general bath spectra, as long as the interaction is bilinear, and the spectra 
ensure relaxation. 

The internal energy reads 



ZTT 



27r 



Twji , / l3Tiuji 



27r 



27r 



27r 



27r 



27r ^ V 27r ) 



(3.45) 



3. The effective temperatures 

We shall now study two objects, — a(x^) and Tp — (p^)/2to, that would in classical equilibrium be equal to T 
and which we shall interpret below as effective temperatures. As in the classical situation, it holds that ||3^ 

= a{x^) = 2a^ (3.46) 
oa 



We find 

ha 



..^ (^i-r)^(^) (c.2~r)^(^) ^ (^3-r)^(^) ^^^^^ 

nm I (cj2 - t^i)(t^3 - t^i) {uji - ll!2){lu3 - LU2) (cji - cj3)(w2 - t^s) 



Likewise, 



ip^) dF, n^v \ -iV'(^) -2^(^) c.3^(^) ^ ^ ^ 

m dm mn (tJ2 - Wi)(w3 - wi) (tJi - W2)(w3 - ^^2) {^^i - i^3)[^^2 - ^^3) 

To find the Gibbsian values for 7-^0 one has to notice that 

^(l^)_^(_l^^)=^T + ^n coth ^ (3.49) 

and this yields the standard weak-coupling result known from all the books, 

U = \t, + = T, = Tp^^ coth ^ = hcoo + (7 -> 0). (3.50) 

2 2^ ^ 2 2 2 e^'^^o _ 1 ' ^ ' ^ ^ 



a. Thermodynamics and effective temperatures at high T 

Using that for small z 



Inr(z) = -In^ -7_EZ + — (3.51) 
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where 7^ = 0.5772156 is Euler's constant, one gets the free energy 



Fp = Tin l3huJo + [r^ — lljI — Lul ~ oj^] w Tin ptiLUo ■ 



24:171 



(3.52) 



where luq = y^a/m and T has been taken large in the second identity. The internal energy and entropy become 



Up = T 



12m 



Sn — In 



T 



6m 



(3.53) 



From eqs. (3.47), (3.4S) we obtain at large T 



12m 



720m2 



(3.54) 



and 



^ _ ^ ^ (3h'ia + ^T) f3^n^ , p^h'l^T^m - (a + 7^^] 



12m 



m 



720m2 



(3.55) 



b. Thermodynamics and effective temperatures at low T 



Further results can be obtained with the improved Stirling formula 

-)\nz - z + - ln(27r) + — 

2' 2 ^ ' 12z 360z 

One gets for arbitrary and for large F 



Inr(z) = (z - -)lnz - z + - ln(27r) + — - (3.56) 



F„ = — [F In r - y iUk \nuJk]~—T^ 

k=l 

[In ^ + 1] - [(1 + w) ln(l + w) + (1 - u;) ln(l - w)] -^T^, (3.57) 



27rm 7 47rm 6ha 

f/,.A[FlnF-^.,ln..] + ^r^ 



^ [In +l]---^[(l + u.)ln(l + w) + (l-w)ln(l-u;)] + — (3.58) 

27rm^ 7 ' Anm^^ ' ^ ' ^ ' ^ " Qfia ^ ' 

Sp = /3(C/p ~Fp)^^ T, (3.59) 
Sna 

Notice that Sp, the shift in total von Neumann entropy due to the presence of the Brownian particle, differs strongly 
from the von Neumann entropy of the particle itself, which remains finite at T = 0, as we shall show in Eqs. (4.34) 
and (4.35). The non-additivity of entropies encountered here is a deep aspect of quantum physics, where a subsystem 
can have a larger von Neumann entropy than the full system. 



Using (3.33) one finds at low temperatures and arbitrary am/7^, w = \/l — Aam/^"^: 

Tp = ^ln ^ + T-^[(l - ln(l - w) - (1 + ln(l + w)] + 0{T^), (3.60) 
■Km 7 47rm w 

T. = ^ilni±^ + 5^T^ + 0(T^), (3.61) 
777 w 1 ~ w Ana 



The above expressions simplify in the limit of strong damping, see Eq. (3.31): 

Fp{T) = 7, (In + 1) + — In ^ ^ , 

ivrm 7 27r7 am ona 



?i7 /, Fm , ha , 7r7 
UpiT) = ir^i\n— + l) + —\n^ + -^T\ 3.62 
27rm 7 27r7 am ona 
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and 



h"/ Tm ha , A, , , 

r=^ln + — +0{T'^), 3.63 

TTTO 7 

T^^—\n^ + ^T^ + 0{T% (3.64) 
7r7 am ina 



4. THERMODYNAMIC ASPECTS OF ADIABATIC CHANGES 

1. Generalized thermodynamic formulation 

We now make clear that the relation with standard thermodynamics can be continued much further by introducing 
the two effective temperatures 

Tp = ^, T,=a(a;2), (4.1) 
m 

One reason to do this is that the stationary state for the harmonic potential has a quasi-Gibbsian expression for the 
Wigner function 

W{p,x) = W,{pW.{x) = -== (4.2) 

with /C(p) = p^ /2m the kinetic energy and V{x) = ^ax^ the potential energy. This expression is quasi-Gibbsian, since 

there occur two different temperature-like variables. (Notice that the normalization is J dxdpW — 1). There occur 
the Boltzmann entropies of momenta and coordinate, 

5p - - 1 dpW{p) HW{p) Vn]^\\n'^ + i (4.3) 

S., = - j AxW{x)HW{x)^] = iln^ + i (4.4) 

(in Sp and Sx we skipped terms In 27r). The complete "Boltzmann" entropy is 

f 1 mT T 

Sb^Sp + S^ = - dpdxW{p, x) \n[Wip, x)h] = - In — + 1 (4.5) 
J 'i h a 

a. Internal energy and interaction energy 
The energy of the central particle reads 



^ + W) = '-Tp + k 
2m 2 ^ ' 2^2 

For a discussion of this identification in systems with out a self-interaction term, see section 

The interaction energy, i.e. the energy of the cloud of bath modes that surround the particle, is defined as 

Uint = C/tot - Ub{j ^0)-U ^Up-U ^Up-^Tp~^Tx (4.7) 



At high temperatures one gets from Eqs. ( p.53| ) and ( p.54[ ) 



Since the energy of the cloud is involves fi,, the non-triviality of the cloud is a quantum effect. 
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At low temperatures one gets the internal energy 

U = -^In^^ - -^[(1 + w)ln(l + w) + (1 - ln(l - w)] + (4.9) 
ZTrm 7 27rm Qna 



For large damping this reduces to 



Ti'y , mV ha „ 7^ , 777 

U=—^lii + —[In-^ + 1] + —^T^. 4.10 

7rm 7 27r7 am b/ia 



The interaction energy of the cloud is now independent of w, 



provided that F is large. 



(7i„t = :/^+0(r4), (4.11) 
27rm 



b. Generalized free energy and the first and second law 



The definition of the effective temperatures admits a clear thermodynamical interpretation. For studying the role 
of an adiabatically slow variation of an arbitrary parameter, such as a or to, that we shortly denote by a, the free 
energy F is defined as 

F^~Tp\nZp-T,\nZ,, (4.12) 
The definitions Zp = J dpc:xjp[—IC{p)/Tp\, Z^ = /da;exp[— V(a;)/r^] bring 

F = -^Tp In TOTp- ir. In ^ (4.13) 

For considering changes in system parameters one needs 

d[-Tp\nZp] = - In ZpdTp - ^dZp = ^ ( dpe-^-^^P'^d]C{p) - SpdTp, 

d[-T,lnZ,] = -InZ.dT, - ^dZ, = ^ f dxe-^^^^^'^dVix) - S^dT, (4.14) 



Eq. ([4.12D then yields 

d^ ^ -SxdT^ - SpdTp + dW,cv. (4.15) 



with, in agreement with the derivations (4.21) and (1.24) below, the work added to the system 



dm da 

dWrcv ^-Tp— + T,— 4.16 
2m 2a 

These relations are valid in spite of the fact that both and Tp are functions of a and m. 



Because of Eqs. (|4.6|), (4.3) and (4.4) the definition (4.12) is compatible with the standard identification |23U2 



F -TpSp-TA (4.17) 

that one would write down immediately for a two-temperature system. From this relation one will indeed reproduce 
the standard formulation for the first law for situations with two temperatures, 

dC/=dQ,ev+d>V,.cv, (4.18) 

dQrcv = Tpd5p + r,d5„ (4.19) 

where dQrev is the heat reversibly added to the particle. A detailed discussion concerning the general definitions of 
the work and heat is given below, in section [7. 1| . 
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The generalized thermodynamical relations (4.14-4.18) are in close analogy with those proposed recently for nonequi- 
librium glassy systems p3|-p5 |. An alogous to that situation, F pertains to the particle alone, and, except at high T, 



it differs from the Fp of Eq. (3.44) in previous section, which relates to the whole system, to be more precise, to the 



particle and the cloud of bath modes around it. 

Let us recall that Fp satisfies Gibbsian thermodynamics, while F does not. There are many physical systems, such 
as a Josephson junction strongly coupled to the electromagnetic field, where the natural object to study is nevertheless 
F, since it relates for that case to properties of the junction only. 

It is standard wisdom that energy is dispersed if the variations are non-adiabatic changes. This is confirmed by Eq. 



(7.43), which holds provided the whole time-domain where m and a vary is accounted for. This leads to the general 
result 



dQ < TpdSp + TxdSx 

that is also known from the study of glasses and, more generally, from two-temperature systems. 



(4.20) 



2. Violation of the Clausius inequality 

a. The Clausius inequality at small T. 

Let US now consider two concrete examples, and study the Clausius inequalitydQ < TdS^N, which is one of possible 
formulations of the second law. 

For a very slow variation of the spring constant a one gets 

dxdpW(p,x) — da^ dxW(x)-x'^da^T^ — , (4.21) 
oa J 2 2a 

in agreement with Eq. ( |4.16| ). The first law implies for the heat added adiabatically to the particle at low T 

dQ_=dC/-dW,.ev = (?^ + ^--)^ = -:^T2da +0{TUa) (4.22) 

oa oa a 2 Sna-^ 



It is seen that dQrcv = at T = for all 7. Using Eqs. (4.1, 4.1) we derive for large 7 and very large T, 
da 



'^^^ 2a 



\n{j^/am) 72ln(rm/7) ^ti^ a'^ Ini-y^ / am) \n{j^/am)' 



(4.23) 



At T = the Clausius inequality says that no heat can be taken from the bath, at best heat can go from the central 
system (here: the Brownian particle) to the bath. In our situation dQ is of order T^, while TdSyff is of order T. 
Since we only do powercounting in T and both expressions are non-trivial, we may replace here SyN hy Sb- Thus for 
the case da > 0, where an amount of workdWrov > is done on the system, the Clausius relation is violated at low 
but non-zero T. 

In the same way one can consider the variation of the (effective) mass m. Here one has 

dWrev = J dxdpWip, x)^dm = ~ j dpW{p)^dm = -Tp^, (4.24) 
again in concordance with Eq. ( 4.16| ). This implies 



dSrcv = (Tp-^ + T,-^)dm = (-i^ H- + ^)— = -^dm + 0{T^) (4.25) 
Om Om om am m Z zttto^ 

In contrast to the previous case there is a transfer of heat even for if the bath temperature is zero. Thus, violation 
of the Clausius inequality is even stronger in this case, since for dm > one has dQ > 0, even though TdSyN — 
(for T — > 0). This situation with dWicv < corresponds to work performed by the system on the environment. To 
underline that the heat comes from the cloud of bath modes, we notice that the general relations 

dWrcv = dFp, dU = dUp - dC/i„t = dQrev + d Wrev (4.26) 

imply 
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dQ,ev = TdSp - dt/int (4.27) 

For changing m in the T — situation it indeed holds that 

dQrev = -df/int (4.28) 



for all values of w, even when T is not very large. For a change in a it holds that dC/int = 0{T ), but Eq. (4.27) 



nevertheless reproduces Eq. (4.22), because of relation (3.59) 



Let us briefly discuss consequences drawn from the violation of the Clausius inequality in the quantum regime. 
First of all it seen that it occurs in the overall Gibbsian state, so that globally (i.e., when applied to the overall closed 
system) thermodynamics is valid by definition. In particular, since the overall system does not absorb heat during 
any variation of a parameter, and dQ = is consistent with the T —^ case of the Clausius inequality (later we 
will see that this is also the case at finite temperatures, where dQ is still zero). Nevertheless, the local state of the 
particle is not Gibbsian and does allow violations as we have seen. We stress that this violation arises due to quantum 
entanglement, which leads to non-gibbsian effective temperatures for the stationary state of the brownian particle. If 
the effective temperatures for T — would equal to their gibbsian values hujo/2, the state of the particle would be 
pure, which is impossible since it does interact with the bath. 

When later discussing the Thomson's formulation of the second law, we will see that it is perfectly valid for the 
overall Gibbsian state, so that the above violation of the Clausius inequality provides us with an explicit example 
showing that at low temperatures the very equivalence between different formulations of the second law is lost. 

A further aspect of this matter is the squeezing of phase space and entropy, relevant for computing in the quantum 
regime. In a separate paper we have shown that the so-called Landaucr bound for the erasure of one bit of information, 
that arises from the Clausius inequality, is violated in a similar manner pj|]. 



Notice again that the effective temperatures remain finite in the limit T ^ (see Eqs. (3.63, 3.64) and Fig. 4.1), 
and both are larger than the bath temperature T. The fact that they are non-equal is due to a mixed state of the 
particle. Indeed, a quantum system non-weakly interacting with its environment, will be in a mixed state even if the 
whole closed system (the particle and environment together) is in a pure state (e.g., the vacuum state). 




FIG. 4.1. The effective temperatures Tp, Tx versus the bath temperature T for two values of the dimensionless damping. For 
the involved parameters we take the following values: Ti-y / (Airm) = 1, TiV /{2-k) — 100. Left figure: am/^'^ — 80 (underdamping); 
from the top to the bottom: Tp, Tx, T. Right figure: the same, but with am/^'^ — 0.2 (moderate overdamping) . 

The existence of different temperatures T^, Tp and T for the subsystem and bath seems to contradict the zeroth 
law, which states that systems interacting for a long time are in equilibrium, and share common temperature. Notice 
especially that the above difference between temperatures is not a consequence of any metastability and/or incomplete 
equilibration, so that our effective temperatures do not depend on the dynamics of the particle and have somewhat 
more definite status compared with those defined, e.g., for glassy systems p^ , p^ , p6[ . Typical derivations of the zeroth 
law (see e.g. Ref. ^ for one of the most clear examples) essentially use the assumptions that the interaction with the 
bath is very weak, and that the total entropy can be considered as the sum of entropies of the subsystem and the 
bath. Evidently, this last condition is not satisfied in our case, except for the limit 7 — > 0, where and Tp tend to 



their Gibbsian value, see Eq. (3.5C), of the harmonic oscillator coupled very weakly to its bath, the situation treated 
in all textbooks. 

Let us notice that in literature some other ways were proposed to establish effective characteristics for non- 
underdamped Brownian oscillator. In [pd[ it is shown that there is a mapping to the Gibbsian (underdamped) 
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oscillator through the definition of an effective mass and frequency. For the description of Hcs Prokof'ev studies a 
related model with a gap in J{u>) at small uj In this approach he makes a different identification for the effective 
temperature. Surely, the choice which quantity to take as "effective temperature" is to some extent a matter of taste, 
that can only be justified by the induced simplification of the physical results. For our thermodynamic approach 
other definitions of effective temperatures will not be very helpful. Our Tp and T^, however, allow to formulate the 
generalized Clausius inequality and they also occur in the Maxwell-Boltzmann-like form ( |4.2| ) of the Wigner function. 
Last but not least, these effective temperatures enter in the same way as in glasses and other two-tcmpcrature systems, 
such a black holes. 



b. Von Neumann entropy 



In next part we shall discuss the von Neuma nn e ntropy of the central particle. To investigate it one needs the density 
matrix corresponding to the Wigner function (4.2). For the harmonic oscillator this can be worked out explicitly. One 
approach is to introduce an effective mass and an effective frequency |ll|] , and insert these results in the expression 
for the entropy of the effective harmonic oscillator. We found it more insightful to redo the derivation. The standard 
relation 



u u 
{x + -\p\x- -) 



Ap e-'P^'/^Wip, x) 



(4.29) 



connects the density matrix in coordinate representation with the Wigner function. From this relation one gets the 
following formula |3q] 



':\p\x') = 



1 



:exp I 



(x + x'Y (x - x')"^ 



(4.30) 



The physical meaning of Eq. ( 4.30 ) is clear: The diagonal elements {x — x') are distributed at the scale ^ (x^) , while 
the maximally off-diagonal elements {x — —a;'), which characterize coherence, are distributed with the characteristic 
scale Ti/ ^ (p^). 

We have to find eigenfunctions and eigenvectors of this density matrix 



Ax' {x\p\x)Jn{x') =_p„/„(x) 



The solution of this problem uses some tabulated formulas for Hermite polynoms, and results in 

1 ■ 



1 



Pn 



V - 



V — 



V - 



Ux)=cH,,{cx)e"'"-"'^, 

I 2\ \ 1/4 



fi^(x^) 



ApAx 



(4.31) 

(4.32) 
(4.33) 
(4.34) 



where 77„ are Hermite polynomials, and it holds that v > ^ due to the Heisenberg uncertainty relation. The result 
for the von Neumann entropy now reads [Tl|] 



SvN = - ^Pn hip„ = {v + ]r) ln(w + \)-{v-]-) \n[v - i) 



(4.35) 



The first terms in its large i>-expansion are 



SyN = In u + 1 - 



1 



1 



1 



24i;2 320^4 2688i;6 



(4.36) 



From (4.3), (4.4) and (4.5) we notice that the same quantity v governs behavior of the Boltzmann entropy 

l = lnt; + l (4.37) 



1 {x^){p^ 
5. = -ln^^ 
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This appears to coincide with the leading terms ol (4.36). It is known to be larger than the Von Neumann entropy, 
and this is obvious from the sign of the correction terms. 

If some parameter (a or m) is varied, then the derivative of SyN with respect to it reads: 



dSyN = In r dv (4.38) 

In other words, the sign of the change in SvN is determined by the sign of the change in v. This holds as well for 
the change in Sb, so qualitatively they carry the same information, and this already was used above to simplify one 
point of the discussion. 

Let us stress that von Neumann entropy SyN is the unique quantum measure of localization, whereas the entropies 
Sp and Sx characterize localizations of momentum and coordinate separately. The differences between Sb and SyN are 
due to the fact that in quantum theory momentum and coordinate cannot be measured simultaneously; in this sense 
Sp and Sx characterize two different measurement setups. Nevertheless, for the harmonic particle if Sy^ increases 
(decreases), then Sp + Sx increases (decreases) as well. Notice that the real importance of Sp and Sx becomes clear 
when they have to be used to generalize the Clausius inequality. The von Neumann entropy cannot be used for this 
purpose whenever Tx Tp. 



c. Clausius inequality at large T . 

At low T only power counting in T was needed for showing the violation of the Clausius inequality. The precise 
definition of entropy, and the quantitative difference between the Boltzmann entropy and the von Neumann entropy 
were not essential, since then TdS —* anyhow. Here we wish to show that the same violation already happens at 
arbitrarily large temperature. To do this we have to use the von Neumann entropy of the subsystem. 



In this section we consider very large temperatures, T ^ HT. Using (3.54) we find from (4.25) for a change in m 
that 

dQ= [1-^i— ]— dm, ASyN = [l-^--r-^ ^^]—dm, (4.39) 

12m 2m 12m 2m 

So, for a change dm > it is seen that at arbitrarily large temperature the Clausius inequality is violated. The 
relative violation (dQ - rdS'„Ar)/dS is of order h^jT/{mT^). 
For a change in a we find 

dQ ^ 1-^ + ^ - ^'^>+ifhm Tda (4.40) 
1 2a 24m 480m2 J ^ ' 



\ 2a 24m 480m2 J ^ ^ 

These expressions differ at relative order /rrv^ , and the Clausius inequality is violated for da > 0. 

The important conclusion of this section is that the violation of the Clausius inequality already occurs at arbitrarily 
high temperatures. Later we point out that a similar conclusion can be drawn about the violation of the zeroth law 
at large temperatures. 

We stress that the Clausius inequality is violated for any finite coupling, and the violating terms only disappear in 
the weak coupling limit 7^0, or in the classical limit h ^ 0, equivalent to the high-temperature limit. 



d. Clausius inequality for comparing two systems. 

For non-equilibrium systems the question has not been settled whether the von Neumann entropy is the true physical 
entropy. As we are inclined to believe that it is, we have discussed that entropy above. 

Let us, however, now consider cases where there is no doubt: For systems in true Gibbsian equilibrium the proper 
entropy of the subsystem is surely its von Neumann entropy. We can now compare two such equilibrium systems, hav- 
ing slightly different system parameters. In standard thermodynamics such a comparison does not yield a new insight. 
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as the equilibrium state of the system is independent of its history. We should point out that in the thermodynamics 
of glasses it is customary to compare cooling experiments at different but fixed pressures. A related comparison was 
also made for black holes: it could be shown that comparing the situation of a single black hole before and after a small 
amount of matter was added, is analogous to comparing two different black holes with slightly different masses p2[ . 
This universality pointed at a thermodynamic behavior of black holes, and the physical framework could indeed be 
provided by one of us, by drawing an analogy with the thermodynamics of glasses Q. 

For our present case we can compare two equilibrium systems at slightly different temperatures. This has the 
benefit that the work term is absent, thus needing no interpretation, and it implies dQ — dU. Using the fully exact 
expressions for the energy and the von Neumann entropy, it is then straightforward to show that at large T 

dQ - TdS^N = ^^ttP- dT (4.42) 
24m 

The standard expectation that this should vanish is again seen to be invalid, and the Clausius inequality is violated 
for dr > 0. As before, the terms in the right hand side vanish only in the weak coupling limit 7^0, the classical 
limit ?i — > 0, or the infinite temperature limit P —> 0. 



e. On our identification of the energy of the subsystem 



In section p. 1| we have considered two physical situations. In the first case the Hamiltonian contains a self-interaction 
term ~ . For that case the above results on the Clausius inequality apply unambiguously. In the second case there 
is no such self-interaction, but the potential energy \hx^ is split as ^ax^ -I- ^"fVx^ ^ and the first part is counted in 
the energy of the subsystem, while the last part is counted with the interaction energy. Let us now shortly look at 
what happens when this is not done, and U — (H) is considered as energy of the subsystem. At large T one will 
have if = U+^{jr/a)T = T + ^{-fT/a)T. Since the work is not modified by this identification, one will have a shift 
in the change of heat, dQ = dQ + ^'jTldT/a — Tda/a^]. From (4.41) and ( 1.42| ) it is seen that then even at very 
large temperatures dQ — TdSvN will not vanish whenever 7 is non-zero. Thus, when there is no self-interaction our 
identification of Ti. as the Hamiltonian of the subsystem is already mandatory for having a proper classical limit. The 
underlying reason is that the Wigner function has the Maxwell-Boltzmann form exp(— [2^^^ -I- ^ax^J/T), involving a 
and not b = a + "fV. This fixes the entropy, and leaves one consistent choice for the energy. 



5. EXACT DYNAMICAL SOLUTION 



We now consider the situation where our closed system starts at time t = 0~ from a Gibbsian initial distribution. 
It could arise if long before the total system was coupled to a 'superbath', that allowed relaxation to equilibrium, 
after which the connection was cut | |5C[ |. A more realistic situation occurs when the bath has small non-linearities, 
that drive the whole system to its global Gibbsian state. 



1. The case when the initial state is a modified Gibbsian 



We assume that for t < the system is in a Gibbsian state at temperature T with certain parameters a — ao, 
771 — niQ, 7 = 70. At t — these parameters are instantaneously changed to a, m and 7, and the system will relax 
to a steady state. This new setup generalizes previous studies with for times t < particle and bath are uncoupled, 
described by 70 = 0. An important benefit is that in the strong damping limit the present initial state can be close 
to the final state, which is, of course, impossible if 70 = but 7 is large. When making the change 70 ^ 7 at i = 0, a 
amount of work 5(7 — 7o)r(a;^)o has to be supphed to the system. This was truly large in our Letter where we 
took 7o — 0, but r and 7 large. In the present setup we can choose 70 — 7, but ag ^ a, implying that the work need 
not be not large, even when the Debye frequency T is large. 



a. The eigenmodes of the initial state 



In the most general case the Hamiltonian has for t < the parameters ao, 70 and mo- It reads 
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with 



" XI ^ + 2 X V^^^ XiAijXj (5.1) 

i>0 * i-j>0 



I p (0) 

Aoo^ °" ^" ; Af,, = Ao = -^ , A,■=a;f<^^, (5.2) 

Too y/morrii 



where c|°'' is given by Eq. (2.20) with 70 replacing 7 



Let us denote the eigenvalues of A by j/^. From a previous section, Eqs. (3.10,3.12), we infer that the eigenfrequencies 
are shifted, 

Uk =ujk- -(fyoii^k) A (5.4) 

TT 

where (po, satisfying < (^0 < ti": is given by 

7or^J^ 

<^o(j/) = arctan 2 , T^2^ 1 (^■^) 



In later sections we only need that for large T 



sin'^o(i^) = 77- _ , .o.,.i/n (5.6) 



The eigenvectors are 



(0) 

eS = «fc, e.^= (5.7) 



with normalization factor 



— = 1 + V ^^3^ (5.8) 

"fe ^TOoTO,(a;/- 1^2)2 

The following normalization conditions hold: 

Y^e^e)^5,,, Y^e'^ei^Ski (5.9) 

k i>0 

For small A one may use 

°° 1 1 

[(i- fc)7r + 0o(t^fe)]2 ^ sin^ M^^k) ^^'^^^ 

to find 



/2ATOo(r2+C^2) /r2+w2 2Aw,siT.^^j 



where (/)q = (j)o{ujk)- In the zero coupling limit 70 ^ one has (j)Q 0, so that i^k ojk and indeed 6ik- The 

latter setup occurs in the standard treatments where bath and subsystem are initially uncoupled. 
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b. The noise 



We consider Gibbsian equilibrium for t < 0. Let us now introduce the creation and annihilation operators 6]., bk, 
satisfying [6fc,6j] = 4; by 



I p = E 



?iA(r2 + vl) sin</,^-moi^fc(z5|e*^'=*-z6ie-''^-*) (5.13) 
7r7or^j/fc 



They indeed satisfy = ihSij due to the normalization condition, as well as [x^p] — ih, [x,pi] = [xi,p\ — 0. For 

t < the Hamiltonian then reads 



k k 

In the Gibbsian state that describes our closed system for t < 0, the density matrix is 

1 



-I3H 



It has the Bose occupation numbers 

{blbk + bkbl) = 1 + 



2 1 

— ^ = coth(-/3nz.fc) 



Combining ( 2.16 ) and ( 5.12 ) we now have for the noise 



i>l k 



(bl + bk) cosLJit + i — ibl ~ bk) sinujit 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



To carry our the z-sum, we have to evaluate 



CiC'l cos UJit 



2jT^A 2Asin0g 



i>i 



^(r^ + "k) 



) cos 



(5.18) 



Gradstein & Rhyzhik present on page 40, Eq. 1.445.2 the equality 

cos kx 1 TT cosh(7ra — ax) 



E' 

fc=i 



2a2 2a sinh'; 



(5.19) 



According to them it holds for Q < x < 27r, but it actually only holds for Q < x < -k, while further it is symmetric and 
periodic. We have to apply this with x — tA, which is surely between and tt, for the cases a = iv j A and a = T/A 



ECi^l cos UJit 
/mn ■ 



H^' + '^l) 



' 27r2A 

A^' + ^l) 



^cosh(7rr/A — Fi) VkCOs{nvk/ A — v^t) 



sinh(7rF/A) 



sin(7rz^fc/A) 



[F sin (l)^ e ^* + Vk cos(0o + ^kt)] 



(5.20) 
(5.21) 



In the last step we have neglected terms of order exp(— 2F/A), which are extremely small. The primitive of this 
relation yields 
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When we insert this in expression (5.18~ 



) for the noise we have the exphcit result 




bk) sini 



-rt 



+ e" 



(5.23) 



The memory of the initial state (i.e. the dependence on Uq, tuq and 70, coded in 0q) is washed out after a time 
T„ — 1/r, apart from a h armless phase factor. Notice that the time-dependencies are exp(ji/fei), as one would have 



expected from E g. (^.12 ). Also note that the exp(— Ft) terms underline the special role of t = 0, just as it does 



elsewhere in Eq. ( 2.18 ). We shall later verify that this passes a consistency check. 



c. The noise correlator 



The noise correlator now decomposes in two parts 

Kis, t) = i ( rj{tUs) + riisUt) ) ^ Kois ~ t) + K,is, t) (5.24) 

The first term is the stationary noise know from the situation where system and bath were initially uncoupled, 

1 f°° 

Ko{s-t)^- dLuKo{uj)cosLj{t- s) (5.25) 
Ti" Jo 

with spectrum 

tanh ^phuj 

It indeed does not involve parameters of the initial state. The connection between properties of the noise and the 
friction kernel is the consequence of quantum fluctuation-dissipation theorem P jT^jll| ] . 

As shown in Ref. [Q, the quantum noise has correlation Ko{t) — — In(rt) > at small times, t <C l/F. At T = 
there occurs for large times the celebrated powerlaw, anti-correlated decay 

Koit)=-^ (5.27) 
This is cut off at times larger than the universal quantum coherence time Tfi = h/T, where 

-12 



(5.28) 



The divergence of this expression at t = shows that a regulator like T is needed. 

Let us briefly explain qualitative reasons for the above structure of the quantum noise. As seen from Eq. (2.16), 
the quantum noise is just a weighted sum of the unperturbed coordinates of the baths oscillators. For T — > the 
unperturbed bath appears in its lowest energy level, and since energy and coordinate do not commute (just because 
coordinate and momentum do not commute), the quantum noise fluctuates even for T — s- 0, and brings a non-trivial 
structure to K{t) in contrast to the classical case, where the noise is just absent for zero temperatures. On the other 



hand, the total intensivity of the quantum noise is zero for T ^ 0: as seen from ( 5.26| ), / dtK{t) ~ 2'yT. For the 
total integral to be zero, the correlator K{t) should change its sign at some intermediate time t. For longer times 
the quantum noise is anticorrelated. The correlator displays a power-law behavior, since the correlation time-scale 
Ti/T is now infinite. A colored noise generated by the low-temperature quantum thermal bath will be the main cause 
of our effects. The classical white noise situation Ko{t) = 2jT6{t) is recovered by taking the high-temperature limit 
(r> ?iF). 

The second term is due to the initial correlation of particle and bath. 
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Ki{s, t) = ifne-^(^+*) + ifi2(s)e-^* + i^i2 ^e"^' (5.29) 

1 roc _ p2 I 2 

= dcji^oC^) o sin^0o(^) (5.30) 



1 /""^ - r 

Ki2{t) = — I duj Kq{uj)[— cos{(j)Q{u) + ujt) — sin{(l)Q{u) + ujt)]&iii(j)Q{uj) (5.31) 
Jo ^ 

where (/)o is defined by (^). 

The standard case of initially uncoupled Brownian particle and bath is recovered for 70 — *■ 0, where (po — *■ 0. 
Then K12 and Kn vanish, making the noise correlator time-translation invariant. In the general case, the initial 
correlation only affects the very short time regime i < r,, — l/F (remember that we assume that T is larger than 
other characteristic frequencies of the damped Brownian particle) . 



d. Variances and covariance arising from the initial state 

The Gibbsian initial state leads to three coupled Gaussian random variables: the random initial conditions, zi = pq 
and Z2 = xq, and the noise z^it) — rj{t). More precisely, when we discretize the time-axis in points i^, the function 
2:3 (t) becomes a set of variables z^^i. Their correlations and cross correlations are 

(p2)=morp(ao,mo,7o), (xg) ^ K{s,t) ^\{r,{s)T^{t) + r,{t)r^{s)) , (5.32) 

{poxo)^0, Si{t)^^{ij{t)po+Pov{t)), S2{t) ^ ^{ij{t)xo + xovit)) (5.33) 



The new terms read 



Si{t) — f di'z/sin(^o(i^) sin((/)o(i^) + vt) coth-/3?iz^ 

t^VTo Jo 2 

(5.34) 

S2(t) — — — — / dvsin(j)oM cosicjyoiiy) + vt) coth-Phv (5.35) 
7^^/7o Jo 2 



2. Exact solution of the Langevin equation 

Next we consider the exact solution of the quantum Langevin equation: 

i = — , p = -ax - j(t)x(0) - [ dt'-f(t-t')x(t') + T](t), (5.36) 
m Jo 

The general solution of Eq. ( 5.36| ) is obtained with help of Laplace transformation. The reader may recall the fol- 
lowing standard relations between functions A{t), B{t) and their Laplace-transforms C{A} — A{s) — dt e~'**A(t), 

C 1^ dt'A{t - t')B{t')^ = i(s)B(s), C |i} = -A(0) + si(s), (5.37) 

One gets 

■msx{s) — mx{{)) ^ p{s) , sp{s) — p{0) — —ax{s) — ^{s)sx{s) + fi{s) (5.38) 

Together this yields 

x(s) = — [mx(0)s + p{0) + 77(5)] /(s) (5.39) 
m 
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where 



ms^ + a + S7(s) 



Thus, the solution of Eq. ( 5.36| ) reads 



xit) = xiO)f{t) + -p{0)f{t) + - / dt'fit - t')r,{t') 



(5.40) 

(5.41) 
(5.42) 



where utQ = y^a/m; f{s) and j{s) are the Laplace transforms of f{t), j{t). Expanding f{s) for large s, one finds that 
/(O) = /(O) = 0, /(O) = 1. Now we turn to our standard case of the Drude-UUersma spectrum: 



For /(s) one has 



fis) 



lit) = 7re 



m(r + s) 



7(s) = 



r + s 



r + s 



r + .s 

(s + r)(ms2+a) + S7r (s + wi)(s + tJ2)(s + cjg) ^ P3(-s) 



(5.43) 



(5.44) 



where Psis) was defined in ( 3.23| ), where also its roots 0^1,2,3 are discussed. 
Likewise one has for the initial Gibbsian states 



fo{s) 



mo(r + s) 



(s + r)(mos2 + ao) + S7or 



One may write 



f^ 



S+iUt' 



(5.45) 



(5.46) 



where, in this connection, loq = lo^, not to be confused with the definition cjo = \/ ajm elsewhere in the work. 
For large F one has 



/i = -h 



m 



h 



(5.47) 



with w defined in (3.32). 



Let us now set, in analogy with (5.13), 



(5.^ 



This implies 



One has from ( |5.4l| ), ( ^.44D and ( |5.46D : /3fe ^ I3{vk) with 



(5.49) 



I3{v) = sin(/)o(j^)E/*(-^* + *'^)e""'* 



7F2i/ 



m(r2 + 1.2) 



sin(/)o(t.)- 



r - cj,: 



(5.50) 



The e terms cancel since, due to (5.44), 
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E 



.n 

Lu, - r 



/(-r) = o 



Next one can check that 



which brings 



^ m(r^ + ;/^) sin ^^^^^^^^ ^ mo(r^ + ;/^) sin0o(:/) 



This leads to the exact resuh 



X e 



sin 00 (i') 



(5.51) 



(5.52) 



(5.53) 



i=l 

/3o(i^) =sin0(zy)e*'^«('')-"^('') 



/3i(i/) = sin 00 (i^) /j 



' — {-LJi + iv) 
70 



7TOo 



m(r + w)(r - Wj) Jam fQ{ii'){LUi + ii') 



(5.54) 



For large times only the first term remains, and the initial condition only enters through its phase factor exp{i(j)o—i4>), 
which has no physical efi^ect, thus showing that the central particle relaxes to its equilibrium state independent of its 
initial condition. The expression for (3i can be simplified by writing it as the ratio of two polynomials, and using the 
fact that uji is a zero of l//(— w), allowing to eliminate the uif term of the numerator. This brings 



f3i{iy) = sin 00(1^) 



a - ao + (mo - m)i^^ f. 



m 



LOi + iv 



which is still exact. Using ( |5.52|) we can also express the result as 



m(r^ + u'^) 



f{iv)e''" + [a - ao - (m - ]h{iy) ^ — ^e"-'* 



i=l 



(5.55) 



(5.56) 



It is trivial to check that, when no change is made at t = (70 = 7, mo = and ao = a), the result = 
sin0 exp(wi) extends the negative time behavior (5.13) to all positive times, even though the noise and the damping 
payed a special (but in that case unphysical) role to i = 0. 

In the rest of this work we shall be mainly interested in the situation 70 = 7, mo = m while oq is different from 
but close to a. One gets in the regime of large T and l/F to linear order in a — ao 



I3{v) = sin0(i/) (,"^^M-'^<P{'^)+^'^t 
where w is defined in ( 3.32| ) and loi^2 in (3.31). 



bJi -\- iv UI2 + iv 



(5.57) 



6. ENERGY OSCILLATION AND NEGATIVE ENTROPY PRODUCTION 

We consider the dynamical evolution of a system initially in equilibrium characterized by a spring constant ao, 
which at t = is instantaneously changed to ai = a. These parameters are connected as 

ao = (1 - Q!o)a (6.1) 

We shall assume that |q!o| ^ 1. We also assume a large Debye frequency F, but this does not lead to principal changes. 
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1. Non-monotonous relaxation of the energy at low T 



Let us now consider how the system relaxes to its steady state. From (5.48) and (5.16) one has 



di^ 



;i(r 



(6.2) 



The infinite time values, discussed already in Eqs. (3.47 3.48| ), can be checked from these expressions. For the 
evolution from the initial state to these values, we shall consider times t ^ l/F, and we can just take the limit F — + cx3 
since no divergencies occur, except for the InF term of the static part (p^) at T = 0. Inserting (5.57) we get to linear 
order in ao 



2 2 27r7 2m 



m 

U{t) 



(p!> 

2m 
2 ^ 



2 P 



1. 



-T 



ha ^ lit 
27r7 2m 
fia it 



with the relaxation functions of coordinate, momentum and energy 



2(1 



Cpir) - - 
w 



Ce{t) = 



w 

oo 



dy y coth(f3h'yy / 4m) 
_^ [il + w)^+y^][{l-w)^+y^ 



^-{l-w+iy)T 



1 — W 



ly 



dy y coth{phjy / Am) 



[(l+l«)2+y2][(l_ 

dy y cot]i{ph^y /Am) 
[(1 — iyY — w"^] 



iy{l ~ w)e-(i-'"+*2^)^ 



1 + w + iy 
iy(l + w)e-(i+"'+*?^)^ 



1 — w + iy 



1 



ly 



(1 - w)e-(i-"'+*?')^ (1 + w)e-(i+'"+''2')^ 



(1 



iyf 



[l+w + iyY 



(6.3) 
(6.4) 
(6.5) 

(6.6) 
(6.7) 
(6.8) 



Of course, one just has Ce = Cjy + Cx- The integration variable is y = Imvj^. The appearance of the dimensionless 
timescale t — ^t/2m is natural, since in the underdamped regime, where w is imaginary, the damping time is just 
Td = 2m/7, see Eq. (3.38). In the overdamped regime the timescales and from Eq. ( |3.35| ) are coded in the terms 
(1 ^ w)t in the exponentials, respectively. In particular, for strong overdamping one has t/r^ = (1 — w)"ft/2m — > at/7. 



a. Classical regime 



At large T the tanh linearizes. One can do contour integration to find for overdamping, i.e. for e = am/j'^ < 4 



and w = y/1 — 4e > 0, the exact results 



ha 



cosh(i 



■ e (1 — w ) 



ha 



sinh(u'r) 



sinh(wr) 



w 



na 



cosh(wr) + 



sinh('i«T) 



+ {l-w^) 



sinhiwr) 



(6.9) 
(6.10) 
(6.11) 



For overdamping {0 < w < 1) these functions are strictly positive. For underdamping one has to replace w — > iw, 

implying coshur — > costBr and sinh(u'r)/u' — s- sm{wT)/w. Then and Cp get zeroes, but remain non- negative, 
while Ce remains strictly positive. 

For the relaxation of the energy this implies 



U = T+^aoT 



cosh(u'T) + 



sinh(u'r) 



w 



+ (1 - w^) 



sinh(wr) 



-2t 



(6.12) 



For strong overdamping (e — > 0, u; « 1 — 2e) this becomes a simple exponential decay. 
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In case of underdamping, e > 7, one has w = iw with w = y/Ae — 1 . This yields by analytic continuation 



(6.13) 



U = T+-aoT 



cos(wt) + 



sin(wT) 



w 



1 + 



sin^(wT) > e 



-2r 



(6.14) 



The term multiplying ao is an oscillating function, and it is strictly positive. Its derivative, 

^aoT 1 + ^2 



U = -- 



m 



sin (wt) e 



-2t 



(6.15) 



has zeros but does not change sign. Physically this means: depending on the sign of aQ, both the kinetic and the 
potential energy oscillate either above or below their final value, and the total energy flow is unidirected; it goes 
towards the bath when a> ao, i.e. ao > and from the bath to the particle in the opposite case. 
In the limit of weak damping {w ^ 1) one gets 



sin(2a;oi) 



U = T+ ^aoT <{ 1 + — v--""^ I e-7t/n 



f7 = _2^sin2(^oi)e-^*/™ 



m 



(6.16) 



Notice that the small but oscillating term in U has become of leading order for U. 




FIG. 6.2. The underdamped situation. Cx as a, function of the rescaled dimensionless time z = t6 = t/\w\, normalized to 
unity at = 0, for large T and S = 0.5 



C 



1 




FIG. 6.3. The underdamped situation. Ca, as a function of the rescaled dimensionless time z = r5 = t/\w\, normalized to 
unity at 2; = 0, for large T and 5 = 0.1 
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6. Weak damping regime 



In the weak damping limit it holds that 



2-v/om 
w = —iw ~ ^ 1 



(6.17) 



For performing the integrals, we write Cx{t) first as an integral from to infinity, and make the shift y ^ w ^ u, 
yielding up to order l/w^ 



°° coth(6 + / it cos ur + sin mt cos (2lu at + ut) 

au- ' 



+ 1 



1 + 



2w 



where a correction factor 1 + u/w in numerator and denominator have canceled, and we denoted 



(6.18) 



b = -Phujo 



Evaluating this to leading order in l/w we get 



C.(r) 



+ coth -/3?iwo cos2tJoi I e 



(6.19) 



(6.20) 



For larg e T t his agrees with Eq. (|6.9| ), of which the last factor now becomes cos^ uiqI. For T — Q and r = it agrees 
with Eq. (3.27) below. Likewise 



Cp(r) = 



7r7 



Ay/am \smh^^phuJo 



coth -/3?itJo cos 2ti>o< I 6 



(6.21) 



As in the classical limit, — C"^ + Cp picks up a contribution of Ce that is subleading but oscillating; it is most 
easily obtained by evaluating in a manner similar to (6.18), 



7r7 



\ . ,2 T ' coth-/3;itJoCOs2a;Qi e 

am \ smh ^Phujf) 2 / 



When inserting Cp + Cx in Eq. (6.5) we have for the leading decay of the energy 



1 



1 



U = -huJo coth -Phuo + -aoT— 



sinh^ ^Phcjo 



(6.22) 



(6.23) 



showing that to leading order in 7 in the weak damping limit the energy does not oscillate, but monotonously leaks 
into the bath (when ao > 0) or is taken from the bath (when ao < 0). At l ow temperature this happens with an 
exponentially small rate. But the rate of energy transfer, determined by ( |6.22 ), 



U 



2m I sinh^ i/3?ia;o tanhi/3?iwo 



cos 2ti;oi e 



(6.24) 



is an oscillating function that changes sign in each period whenever T is not infinite. Thus the rate of energy transfer 
is not uni-directed except for the classical limit. When averaged over one period, the cosine is subleading and a 
unidirected flow emerges. 



c. Quantum regime for non-weak damping 
At T = one has y coth(/3?i7j//4m) — \y\. For time r = one finds by direct integration 

Cx = — - „ 3 In , 6.25 
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and we define the short hand 



It further holds that 



AH-T^ln^^ (6.26) 
^ ' 2w 1-w ^ ' 



Cp(0) = -C',(0), Cb(0) = (6.27) 
These results can be verified using the relations 

a{x\t = 0+)) - T,{a) = -TM = -ao«'%^ = ao—C^O) 

ao da 7r7 

-{p'it - 0+)> - Tp(a) = -aoa^ = «o — Cp(0) (6.28) 
TO da TTj 

So after the instantaneous change of the spring constant {t ^ 1/r), the deviation of the potential energy from its 
final value is, to leading order in ao, just opposite to the one of the kinetic energy. Consequently, the particle has 
already the proper energy, but this will not remain so; for ap > first a flow from the bath will occur and then a 
reversed flow, after which the equilibrium will be reached by a second energy flow from the bath to the particle. 
At large times (r ^ 1) and still T = one gets in case of overdamping {w > 0) 

1 -rr , sinhwT, 

Cx (t) = - ^ ^ e [cosh WT 



2e^T^ w 
1 ^ sinhwr 



Cp(r) 

Ce{t) = --^ e-^ [cosh u.r + ] (6.29) 



and for underdamping (w > 0) 



„ , , 1 _ sinwT, 

Cx (t) = - e cos WT H : — 

2e"^r^ w 

1 _^ siniDr 



Cp{r) 

Ce (t) = ~ -i-^ e-^ [cos ^Dr + ^^^] (6.30) 



The latter expressions all exhibit an infinity of oscillations around C = 0. For overdamping one has Cx(0) > 0, 
while it has a negative tail; consequently there remains one oscillation even in the limit of strong damping. In that 
limit (7 large), one may set cr = (1 — w)t = at/^. For large, but fixed a one gets Cx = Ce = — 2 exp(— (t)/(t^, 
Cp = — 8£exp(— cr)/cr^. 



d. Strong damping at low T 

Let us now write Cx (t) as 

Cxit) = fxir,w)+fxiT,-w), (6.31) 

f - ^iin^.-Ci-"')^ Z""" ^ . . y coth( by/Vl-w^)[{l-w)cosyT-ysmyT] 

w J_^ [(1 + u;)^ +2/^J [(1 - w)^ +r 

where b was defined in Eq. ( |6.19| ). 

We investigate in some greater detail two particular cases: w 1 (strong overdamping) and w = (the border 
between overdamping and underdamping). For the first case one changes the integrating variable y y/{l — w) and 
arrives at: 

rf N f to \ ff\ --y A ycoth(6?/v^)[cosycr - ysinyo-] 

fx{T,w) = fi{2eT), /i(cr)=e / dy , (6.33) 

J —OQ [ ^ ^ y \ 

fx{T,-w) = -ee J^^dy [l + y2]2 • (^'^^^ 



29 



Recall that w = \/\ — 4e, and in the limit e ^ one has w = 1 — 2e — > 1. It is seen that in this limit fx{T^ —w) is 
small compared to fx{T,w) due to an extra factors e and, above all, a quickly decaying exponential e"^"^ — e^"l^ . 
Thus, we will omit jx{j-, —w). Then one has a scaling form: 



CEiT,e)=CxiT,e) = fi{2eT) 



hi-)- 



(6.35) 



Notice also that for this function small and large temperatures are determined by the dimensionless ratio: by^ ~ 
■^(3ha/j = ^h/{TxT). If this parameter is small (which is always achieved for large temperature and also for fixed 
temperature and large damping), then we go back to the situation of eq. (6.13), 

ttT 



/i(2£r) 



-4er 



In the zero-temperature limit one takes t/coth(6y\/e) = \y\ and gets: 

y [ cos ya — y sin ya ] 



/i(a) = 2e- 



dy ■ 



[1 + y^r 



(6.36) 



(6.37) 



This function can be exactly expressed through Meijer functions, but we will not write down this representation 
explicitly, since it is useful only for numerical computations. Notice that (6.47) can be once more checked with help 
of ( |6.37 ). The behavior of /i(it) for different temperatures is presented in Fig. |6.4| . 




FIG. 6.4. The case of strong damping. fi{t) as function of dimensionless time t, normalized to unity ai t — 0, for d ifferent 



values of the dimensionless temperature 9 — l/(6y^) = 2T'y / (ha). Upper curve: 9 ^ oo (the expression given by (3.36) taking 



into the normalization). Middle curve: 9 = 1 (see Eq. (6.33) ). Lower curve: T = 9 = (the expression given by (3.37) ). In 



the latter case there is still an oscillation, despite the strong damping. 

It starts with /(O) = 1, becomes negative at tq — 0.407211889989, goes trough a minimum, and finally bends up to 
0~ for T — > oo. The minimum is characterized by 

fjmin = 0.87908730804, /i(crmi„) = -0.0918980496, ca = /"(crmi„) = 0.404842 (6.38) 

In this limit Cp has an interesting behavior. We discussed already that Cp(0) = —Cx{0)- For small cr, Cp quickly 
grows, goes trough zero, and then becomes of order e, starting as ehul/a for small, but not too small a. For finite a 
one thus has Ce ~ Cx, implying that now the total energy makes one oscillation, despite the strong damping. 

For Qfo > it says that, after initially energy has been put on the particle by the change of ao — > a, this energy 
leaks away into the bath. However, at intermediate times more leaks away than in the final state, so a part has to 
come back at moderately late times. This non-monotonous behavior ("bouncing") is familiar of the noise correlator, 
which is anti-correlated at large times in the quantum regime. 

Let us now turn to the behavior of Cx for w = 0: 

a.(r)=8re- ^ ^"^^^^ ^ 7 ^ U 8 /°° jy V ^"^^^ ~ v') ^^^V- - V- 1 . (5.39) 
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FIG. 6.5. The case with w = (the border between overdamping and underdamping) . Cx{t) as a function of dimensionless 
time t, normahzed to unity at t = 0, for different values of the dimensionless temperature = 1/h — T\/2m/ {ti^/a). Upper 
curve: 9 = 1; middle curve: 9 = 0.2; lower curve: T = 6 = Q. In the two latter cases there is an oscillation. 



For Cp{t) one has: 

Cp{t) = fp{T, w) + /p(t, -w), 



Ay 



y coth( b y I \/\^-w^ )[y'^ cosyr + y{l — w) sinyr] 



W " [(l+W;)2+Z/2][(l_y;)2+y2]2 

The behavior of Cp{t) in the overdamped situation can be studied along exactly the same lines as for Cx{t)^ 

y coth( b j/Y^ ) [ y2 cos(2 e ry) + y sin(2 e ry) ] 



Cp{t)= 2ee 



-2e- 



dy 





oo 



Ay 



[l + y2]2(i + e2y2) 
coth( b y/ Ve ) [ V cos(2 ry) + sin(2 ry) ] 



2 12 



(6.40) 
(6.41) 



(6.42) 



Due to the additional factor e, this is smaller than Cx{t) for r > 1, and this justifies Eq. ( 6.35| ). Nevertheless, in the 



qualitative level Cp{t) displays nearly the same behavior as Cx{t)- This is demonstrated by Fig. 6.6 




FIG. 6.6. The overdamped situation: e = 0.1 (w = 0.7745). Cp(t) as a function of time, for different values of the 
dimensionless temperature 9 = l/{b^/e) = IT-y /(Tia). Upper curve: 9 = 33; lower curve: 6 = 0. 
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FIG. 6.7. The zero-temperature behavior of the C-functions versus the dimensionless time. Left figure: w = 0.9 (strong 
overdamping) , 1: Cx, 2: Ce = Cp + Cx, 3: Cp. Right figure: the same but with w = O.I (weak overdamping) . 



Let us now investigate properties of Cx(t) in the underdamped limit , wher e w = iw = iy/Ae — 1, and w is real. In 
the additional weak damping limit one has zz) ~ I/7 ^ oo. Using ( |6.32 , |6.33 ) one gets: 



; e im 



dy 



ycoth (^—^^=^ [{1/w — i) cos {yrw) — ysin [yrw) 



212 



Im 



{l/w- 



1 1 e 



dy 



V^°^^ ( vffi^+i ) [2/^ cos (yrw) + {l/w - i)y sin (yrw)] 

[{l/w + i)2 + y2] ^^l/^ _ j)2 y2]2 



(6.43) 



(6.44) 



The behavior of these functions, as well as CE{t) — Cx{t) + Cp{t), is depicted in Figs. 3.8, 6.9 for T = 0. It is seen 
that for the initial time of order l/w, Cx and Cp oscillate with the amplitude higher for larger w. 




FIG. 6.8. The underdamped situation. 1: Ca; as a function of the rescaled dimensionless time t = r/\w\, for T = and 
w = 2\ 2: Cp as a function of t; 3: Ce as a function of t both for the same values of the parameters. 
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Cp , 3 \ 



FIG. 6.9. The underdamped situation. 1: Cx , 2: 
T — and w — 10. As compared with Fig. (6.8), the amplitude of oscillations is much larger 
are quite close to each other. 



C_B, as functions of the rescaled dimensionless time t = t/|ui|, for 

It is seen also that Cx and Cr, 



e. Moments of the relaxation functions at low temperature 

We can also determine the integral 

h drC.(r) = 8(1 - u; ) ^ ^ ^,^3^^^ _ ^ ^^,^3 [(1 - ^ ) - 2, (1 + ^ ) ^ 3, ] (6.45) 

At T = it can be simply checked it vanishes at w = and near w = \: 

/ dTC,(r) =8[i?(l,5)-2B(2,4)-3B(3,3)] =0; / dTC,(r) - (1 - / dy\y\j}- ^— 4il = « 

(6.46) 

where B{w, z) — T{wyr{z)/T{w + z) is the Beta-function. It actually holds that for all w that 

/>oo 

/ AtC^{t,T = Q) = Q (6.47) 

This surprising zero-temperature outcome will have important consequences when it comes to work extraction. It is 
similar to dtK{t) — 2jT for T ^ 0, where K{t) is the autocorrelation function of the quantum noise. For 
small T one gets 



It can also be verified that at zero T 

1 ^(1) 



f 1 

Cf^ = - dr rC, (r, T = 0) = — C^^ 

^(1) _ (3 + u;^)(3w^ - 1) + 3(1 - w'^fXjw) 



(6.49) 



where A is defined in Eq. (|6.26D , and 

poo o 1 

These coefficients are positive for all w. The minus signs in the integrals arise because the negative tail of Cx{t, T = 0) 
gets a larger weight than its positive center. These results follow from the Laplace-transform 

(2m, T = 0) - / dr (r, T = 0) e-^-- 
Jo 

1 l-w^ r(l + 2M-w)ln(l + 2u- w) (1 + 2u + w)ln(l + 2m + w) 



2uw^ 8^2(1 + u) 



{u — w)w{l + u ~ w) {u + w)w{l + u + w) 



(6.51) 



2u^(l + uV + 2u(l + u)'w + w^ -w^ , , , 2^2(1 + -u)2 - 2u(l + u)w + + , , 

i -^7^ '- ^ ln(l - w) i -4?^ '- ^ ln(l + w) 

[u -V w)w''\ \ -V u — w) ^ w)w^(l + u + w) 



At T = all even moments of Cp vanish. This implies in particular that the integral of = Cp + C^; vanishes. 
The Laplace transform of reads in that case 



/>oo 

Ce(2u,T = 0)= dTCE{T,T = 0)e- 
Jo 



2UT 



(l + w)(l + 2u + w) l + 2u + w {1 - w){l + 2u- w) 1 + 2u - w 



Au^wi^l + u)(l + u + w) I + w 4:U^w{l + u){l + u — w) 1 — w 
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The Laplace transform Cp(2u, T = 0) follows as Ce(2u, T = 0) - C^{2u, T = 0). 
For later use we introduce the coefficients 



Ce' = - / dTTCEir) = Trr^+ 0{T^) (6.53) 



.(1) _ / A^^n^i^\ - , rnirr2\ 

24e2 



They differ from the Ci"'^'^^ only by the factor Cx''' , which goes to unity for large damping. 



2. Entropy production versus energy dispersion 

To derive the rate of entropy production we first need the Wigner function and its temporal evolution. 

a. Fokker-Planck equation for the Wigner function 



To derive the evolution equation for the Wigner function, we shall write the Langevin eq. ( 2.18 ) in the form 



p = —ax p + rj + Sp (6.54) 

m 



where 6p is small as l/F. Indeed, from this definition and the exact dynamical solution (5.41,5.42) one may derive 



6p{t) = mxag{t) + pag{t) + f dt' g{t - t')r]{t') (6.55) 

Jo 



3 ,2 



where 

m ^ — ' i — uji 

i—l 

is of order l/F for large F. Now recall that for the harmonic situation the Wigner function is given us 



W(j>,x,t)^ J dpodxoW{po,xo,0){S{p{t)-p)S{x{t)-x)), (6.57) 
where the average is taken with respect to the noise, W{p, x, t) and W{pq, xq, 0) are final and initial Wigner functions. 



while p(t), x{t) are the solutions of (2.18) for the corresponding initial conditions, and for a particular realization of 
the Gaussian noise. Eq. ( ^.57 ) is not the most general definition of the Wigner function, but it is exact for harmonic 



systems. 



We now seek a closed equation for the Wigner function (6.57). Differentiating W{yi,y2,t) we get 



dt ^ dyk dyi 



where 



vi — —ax — —p, V2 — —■ (6.59) 

m m 



are the damped Newtonian acceleration and the velocity, respectively. The term 6p is a linear combination of po , xq 
and r]{t'). Due to t h e Gib bsian initial state, these are Gaussian random variables and their cross-correlations were 
given m Eqs. (|5!3^)- (|5!35|) . Let us denote these variables by the vector z — {po,xo,ri[t)}, and their correlations by 



the matrix Mij = (ziZj). One then has for its joint distribution 
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Using the relation 



(6.60) 



(6.61) 



one can perform a partial integrations, which brings a closed equation for W . The final result is that we obtain a 
diffusion-type equation (Fokker-Planck-Kramers-Klein equation) for W itself: 



dw{p, x,t) ^^^^ pdw ^ d 



dt 



m dx dp 



[l.p + ax)w]+[D,{t)-D,{t)\£^W + ^D,{i)^, (6.62) 



where the diffusion coefficients and Dp are instantaneous functions t. (Notice that in ref. g we used the notation 
Dpp = Dp, D = Dx, Dxp = Dx — Dp.) The derivation along this road is somewhat lengthy. A quicker way to derive 
the result is to use the solution of the Fokker-Planck equation, determined by the moments 



(p2) ^ 2mif (i). 



a 



'-V{t) 



(6.63) 



where K{t) = {JC{p)) and V{t) = (V(x)) are the expectation values of kinetic and potential energy, respectively. The 
time-dependent Wigner function thus reads 



Wip,x,t) 



27r V 'ianiKV - m'^V^ 
Inserting this in equation ( 6.62| ), one finally gets 



exp 



aKx'^ — Vpx + Vp'^/m 
4KV-mV^/a 



(6.64) 



Dpit) ^ 2Kit) + - [Kit) + V{t)] 
7 



(6.65) 



777 . . • 777 " 

Dxit) = 2V{t) + -[K{t) + V{t)] + -^V{t) + -Vit) (6.66) 
7 a a 



Let us also define the time-dependent but current-less state 



Wst (p, x,t)^ ^ exp ( - ] (6.67) 

^ ' 27T^mDp{t)Dx{t) ^\ 2mDp{t) 2Dx{t) } ^ ' 



for which indeed the right hand side of the Fokker-Planck equation ( 6.62| ) vanishes, though the left hand side does 



not. This is the locally-stationary distribution. For sufficiently long times, that is when Dx(t) and Dp{t) are changing 
with time slowly enough, Wst becomes a solution of the Fokker-Planck equation. 

b. H-function and entropy production 

The TL-function is defined as the information theoretical distance between the actual Wigner fimction W{x,p,t) 
and the locally-stationary Wigner function Wst{x,p,t): 

f dxdp W{x, p, t) In > 0. (6.68) 

J Wst{x,p,t) 

The 7i-function is non-negative due to the inequality 

WWW 

n> [ dxdp[W{x,p,t)-Wst{x,p,t)]^0. (6.70) 
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Thus, Ti, is equal to zero only for W(a;, p, t) = Wst {x, p, t) , i.e., in the stationary state. Since values of Ti, at intermediate 
times are higher than its final value, it is reasonable to look at its rate of change. In particular, Ti. changes with time 
due to the time-dependence of the reference Wigner distribution Wst, whereas the remaining part of dTi./dt appears 
to be induced solely by the bath (see below). We define the entropy production diS/dt by 



dH 
It 



-7- = I dxdpWst{x,p,t) 



This leads to 



d^ 
dt 



dxdpW{x,p, t) In 



W{x,p,t) 
Wst{x,p,t) ~^ 

W{x,p,t) 



d^S_ 

dt ' 



Wst{x,p,t) 



(6.71) 



(6.72) 



This definition has the following properties: 1) It is equal to zero in the stationary state; 2) It is equal to zero if the 
brownian particle does not couple with the bath; 3) It is non-negative in the classical case, where = Dp = T . The 
last two prope rties are proved below. 
Using ( |6.67 ) and denoting 



R{x,p,t) 



W{x,p,t) 
Wst{x,p,ty 



(6.73) 



one gets: 



m 



dx dp R{x, p, t) Wst {x, p, t) 



dxdp[CW{x,p,t)] lnR{x,p,t) - 
= J dxdpW{x,p,t) \nR{x,p,t)~ J dx dp R{x,p,t)Wst{x,p,t). 
where C is the Fokker-Planck operator of the right hand side of Eq. ( |6.62| ) . Noting that 



CUnR=^C^ R-^Lopit) 



dR 



1 2 



dp 



[D^{t) - Dp{t)] 



dR OR 
dx dp 



and making once more integration by parts one ends up with 



d^ 
dt 



dx dp 



W{x,p,t) 

RHx,p,t) 



iDpit) 



dR{x,p, t) 
dp 



[D^{t) - Dp[t)] 



dR{x,p,t) dR{x,p,t) 



dx 



dp 



(6.74) 



(6.75) 



(6.76) 



Now it is clear that in the classical white-noise limit, where — Dp, the entropy production is non-negative. The 
positivity of diS/dt just means that from the global viewpoint the approach to the stationary state is monotonous. 
In contrast, in the quantum case the positivity of the entropy production is endangered. It is also clear that for a free 
brownian particle (7 = 0) the entropy production is zero. 

Finally we mention that the difference between dSs /dt and the entropy production is just the entropy flux: 



deS 

~dr 



dSs ^ ^ _ ; ^xdpW{x,p,t)\nWstix,p,t). 



dt 



dt 



It takes the value 



deS 
dt 



k{t) , v{t) 



Dp{t) D,{t) 



(6.77) 



(6.78) 



Let us recall that in the relaxation process no work is performed, so a change in energy can only be due to a change 
of heat exchanged with the bath. Therefore the last relation can be written as 



deS _ Qp Qx 



dt 



D„ 



D, 



(6.79) 



where Qp — K and Qx — V are the changes of heat in the momentum and coordinate sector, respectively, while Dp,x 
are the corresponding diffusion coefficients in the Fokker-Planck operator C of Eq. (6.62). Notice that this entropy 
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flow deviates from the standard expression d^S/dt = Q/T = Qp/T + Qx/T, that could not make sense since Q does 
not scale with T at low T. 
The Boltzmann entropy reads 



Sb 

Its rate of change is 



dpdxW{p, X, t) In ^W{p, x, t) = 1 + ^ In 



h^n V a 



dSs '2aKV + 2aKV ~ mVV 



dt 



AaKV - mV^ 



(6.80) 



(6.81) 



The entropy production is the difference between them and appears to be quadratic in the deviation from the 
equilibrium state. To second order in the small parameter aa it becomes 



diS _ h^^aal ( C: 



dt 167r2 



T2 



+^(ci + c:)(77l + 77#) + 



■ CC 



(6.82) 



where C^ .^ denote the dimensionless derivatives dCp^xi^) / dT . 



c. Classical limit 



In the classical limit, where Tx = Tp ^ T, the rate of entropy production this becomes the sum of two squares, 



much alike the energy relaxation function of eq. (6.11) 



dt 7 



coshf 



sinh WT 



sinh WT 



(6.83) 



The total entropy production is thus 



This result holds for all 



dt 



d^ ^ q| 
dt 4 



(6.84) 



In Eq. ( 7.42 ) we derive the general result for the energy dispersion. In the present setup we have a{t) = aoS(t), 
yielding 



An = ^ 

4ttj 



CxiO) 



With help of Eq. ( |6.9[ ) we find 



An _ al 

This just coincides with AjS', explaining that both describe the same physics. 
In the strong damping limit w —f 1 one a simple exponential decay, 

d^ ^ aa| ^_4at/7 
dt 7 

In the weak damping limit, but still at high temperatures, the result oscillates, but is non-negative, 

2 

m 



M^2:^sin2^oie-^-*/™ 
dt 



(6.85) 



(6.86) 



(6.87) 



(6.88) 



Notice the similarity with the rate of energy decay (6.16) 
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d. Weak damping limit at moderate temperature 



In the we ak-da mping limit 7 — > 0, where Tp — — ^?ia;ocoth(i/3?ia;o), the entropy production foUows from ( 3.82 ) 
(S.20) and (6.22) to leading order as 



^ = 2^ I sin^ 2uJot + (cos 2u;ot 
dt 4m ' 



. ^ ^1: — y + - Alp tanh -/3?iwo smAcuot 
smn pnujo n 2 



-27t/ji 



The term Aip comes from the difference Tp — T^, given in (3.4S). It reads 



Alp = ■iPi^^) - ) 



(6.89) 



(6.90) 



At high T the last contribution of (6.89) vanishes (at least, it is of order l/F, which we discard everywhere in this 
work), so Eq. (6.88|) is recovered. 



The term with Atp is responsible for the occurrence of both positive and negative values of the rate of entropy 
production. At low T one has A'tp — > InF/wo, which is moderately large. Therefore, below some specific temperature 
T* ~ hujQ / \ti{T / LiJo) the rate of entropy production can be negative, a surprising result. This finding goes against the 
formulation of the second law in the form of positivity of the entropy production. In our system the negative rates 
are not totally unexpected since oscillatory behavior is also exhibited already in the rate of energy decay ( 6.24 ). 

The integrated entropy production is, to leading order in 7, insensitive to the oscillations. When averaged over one 
period, the cosine and sine are subleading and a positive rate emerges. The full integral reads 



while the energy dispersion is 



A,S- 



An 



1 



1 



sinh^ [BfujjQ J 



(3hujQ 
sinh phcuQ 



(6.91) 



(6.92) 



Both expressions have the same order of magnitude, and coincide at large and small T. 

Consequently, in the Gibbsian limit the rate of entropy production oscillates in case of underdamping, as does the 
rate of internal energy. After averaging over one period the oscillations are wahsed out. This justifies our identifications 
of entropy flux and production. 



e. Entropy production at zero temperature 

Also at zero temperature the entropy production can be negative. Let us consider the case of strong damping, 
where e -C 1 and Tx ~ eTp, implying 

^ - {i^ + ^)CL' + \C'.C':] = ^^/V)[/((-) +e/r(.)] +0(e^) (6-93) 



Now we know that /i has a negative minimum at tTmin- Let us expand, using the numerical constants from Eq. (3.38), 

h{a) = /l(a,„i„) + ic2((7 - CTniin)' (6.94) 

Then we get 

rl l? n(^ 4^ 

-C^(cr - crmi„)(cr - (Tmi„ + e) (6.95) 



djSx _ a(l + e)a l 
dt ~ 47A2 



This is negative for dniin — £ < cr < Umin- The minimum is of order — £ , and the area of the negative part is of order 
£^. Notice that negative value holds over a time window 5a = £, corresponding to At = T.p = m/j. This is much 
less than the free oscillation period tq = \frnfa^ so after averaging over one period it disappears. However, in the 
(strongly) overdamped regime there are no oscillations, so there is no compelling reason to average over one period. 
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For a numerical investigation of the entropy production at T = we will first of all introduce a new parametrization 
for the effective temperatures: 



■Km 



1 

4it; 



A 



l + w 



A 



1 — w 



nm 



1 , l + w 
-In- , 

w I — w 



(6.96) 
(6.97) 
(6.98) 



where A — 2mT/"f is a large dimensionless parameter. Then Eq. (6.82) can be presented in a more convenient form: 

^ - ^ ^ (^^" + ^(c. + + e's'c;) + c^c::) . (6.99) 

As for the functions involved in this expression, one notice: 



dr w 



dy 



ycoth (^by^J\^^ cos[y{l- w)t] 

[ (1 + W)2 + (1 - W)2y2 ] [ 1 4. y2 ] 



dr^ w 



dy 



ycoth (^by^^'j (cos [y(l - w)r ] + ysm[y{l - w)t]) 

[(l+U;)2 + (l-,«)2y2][l+y2] 



dr w 



dy ■ 



y^cothl^by^^) sin[y(l-u;)T] 

' (1 + w)2 + (1 - w)2y2 ] [ 1 _^ y2 ] 



(6.100) 
(6.101) 
(6.102) 



Recall that Cp,^{T) = fp,x{T, w) + /^.^(t, -w). 

The behavior of the entropy production is depicted in Fig. 6.10. It is seen that there is a small region, where the 
curves are negative. For w — 0.1 (weak overdamping) the negative region is yet noticeable, but already for w = 0.7 
(moderate overdamping) this region is almost indistinguishable. This is in agreement with the above analytical 
analysis in the limit e — > 0. 




FIG. 6.10. Rate of entropy production (up to a factor aoi^/ {Id'^e^Ol) ) versus dimensionless time r. Left figure: T = 0, 
w = 0.1. Right figure: magnification of the region around r = 1.75. 



7. WORK AND HEAT 



So far we have discussed the system's relaxation from a non-equilibrium initial state. Since the total system is 
isolated, in this process energy is transferred from the subsystem to the bath, or vice versa. This energy is related 
to the unobservable bath modes, so it is identified as heat. In particular, in no way work was added to or extracted 
from the system, except for the initial moment, where the strength of the central spring was modified. 

We shall now consider the possibility of additional changes in the spring constant and its implications for work 
extracted from the system. 
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1. General definition of work and iieat 



The behavior of a statistical system under interaction with external macroscopic sources is the standard area of 
applications for any thermodynamical theory. As known well, in this setup one neglects the influence of the statistical 
system to the dynamics of the source. Therefore, one can keep the parameters of the system as given functions of 
time, and solve the corresponding equations for the system's dynamics. We start with general remarks about the 
energy budget of any variation. 

Let us consider the change of a system parameter a. It is assumed to be intrinsic, that is to say, to characterize 
the Brownian particle but not the bath or the interaction between the particle and bath. In the situation discussed 
in the body of this work, a can stand for the spring constant a of the harmonic potential, the effective width of the 
anharmonic potential, or the mass m (for electrical circuits and junctions mass is connected with inductivity and can 
be subjected to variations; this also makes sense in systems where m is an effective mass, that can be modified by 
changing other system parameters, such as pressure) 



First one has to identify the Hamiltonian of the subsystem. In Eq. (2.1) we have choosen H — ^p^/m + ^ax^ 



as in absence of the bath. It should be stressed that we did not include the self-coupling ^'yTx^ or (part of) the 
interaction energy in Ti. Our choice is the natural one in the sense that the limits of large Debye frequency F and 
subsequent large damping 7 lead to moderately large values of the energy of system and bath, and not to large terms 
of opposite sign that cancel i n the total energy, as would occur, e.g. if ^"fTx^ were counted for the subsystem. See 



also the discussion in section 4. 2 e 



A change with time of the mean energy is considered due to a variation of a parameter a according to the prescribed 
trajectory a{t): 

au.^l d.dp W(,. .) = / dxdp«(p. .) dinp. X. ,) + / d.dp Wi,, t) d«(p. (7.1) 

where W{p, x) is the Wigner function of the Brownian particle. The last term results from the change in the Hamil- 
tonian, so it is a mechanical, non-statistical object. Following other authors ||,^,^,^,^, we shall associate it with 
the work dW produced by external sources, in close relation with the definition of work in classical mechanics and 
standard thermodynamics. The first term in the right hand side represents the variation due to the statistical redis- 
tribution of the phase space. We shall identify it with the change in heat dQ. Eq. (7.1) can then be written as the 
usual first law: 

AU=m+SW (7.2) 



The work, as defined in Eq. ( 7T ) can be shown to be the change of the total closed system's (the particle plus bath) 
energy due to the variation of the parameter a. First one notices that for closed systems with a unitary evolution 
any change of energy is determined solely by work. This fact is due to the conservation of energy, and can be easily 
illustrated using the von Neumann equation of motion for the density matrix ptot of the total system. Indeed, since 

dptot i 



dt h 



[ptot,'H 



tot 



one has 



AUtot^ d iy{ptotHtot) = tr(pfot dTLtot) + tr{Htot dptot) 

1 

= ti{ptot dHtot) - T dt tr{Htot [ptot^'Htot]) = tr(/3tot dUtot) (7.3) 
n 

due to the cyclic character of the trace. If now Wtot{p, x,pi, Xi, ...) is the Wigner function of the whole system, then 
this implies the identity 

dUtot = dp dxY\dpkdxkWtotip,x,pi,Xi, ...)dntot{p,x,pi,Xi, ...) = dx dpW{p,x,t) dn{p,x), (7.4) 
k 

since we only consider cases where forces are attached to the central particle, so dT-ltot{p,x,pi,xi^ ...) — dTi.{p,x), 
implying that the pk, Xk integrals over the full Wigner function just bring the Wigner function of the subsystem. 
Taking into account that no heat was added to the total system, dQtot = 0, we may conclude that the work dW 
extracted from the subsystem equals the work subtracted from the total system, dWtot = dUtot = dW. 
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We again consider the situation where at t = the system was Gibbsian, with spring constant oq, and it is 
instantaneously changed to a new value a. For achieving this an amount of work Wo, given in (7.13) has to be added 
to the system. For a > uq this is just the work needed to make the spring attached to the particle stiffer. When 
ao < a this work is negative; energy is extracted, since the spring is weakened. For work extraction from the thus 
created non-equilibrium state we shall make additional changes in the spring constant. 



2. Maximally extractable work 

Now the total system consisting of central particle and bath is out of equilibrium, some work can be extracted from 
it. Before the interaction between the bath and the particle has been switched on, the total energy was 

t/tot(0-) = Ub{T) + Up{T, ao) (7.5) 

where Ub{T) = 7r^T^/(6?iA) and UJT) are the initial energies of the unperturbed bath and the perturbation due to 



. . . . p V , 

the Brownian particle, defined in ( ^.17| ) and by Up = d[l3Fp]/dp with Fp taken from (3.44), respectively. After the 
switching of the interaction has been completed, the energy has become 

f/tot(0+) = Wo + C/B(T) + C/p(r,ao) (7.6) 

Let us now consider what is the maximal amount of work that can be extracted from the overall isolated system in 
the considered non-equilibrium state. First of all, we notice that we are interested in the work done due to the non- 
equilibrium character of this state, and not by a work which might be done due a change of the Hamiltonian. Therefore, 
during extraction processes the parameters of the Hamiltonian Titot will be either fixed or vary cyclically, such that 
after the process has been completed, the system has the same Hamiltonian as initially. To determine the maximal 
amount of extracted work we will employ the following formulations of the second law, which are undoubtedly valid 
for the considered thermally isolated system P, |27| , p3[ | . (i) No work can be extracted from a system in its equilibrium 
state. (Let us recall that thermal isolation means that no external supply of heat is given; the allowed transformations 
are variation of parameters by external sources), (ii) The converse is true as well under certain general conditions 
p3| : If no work can be extracted by any means from a system in a given state, then this state is equilibrium. 

As follows from (i) some work can be extracted from non-equilibrium states. In the same way (ii) implies that 
if work has been extracted in all possible ways, the system is left in an equilibrium state at some temperature T^n. 
Since the overall system is thermally isolated, the extracted work is in magnitude equal to the complete change of 
energy (this is the statement of the first law): iVVmaxI = C^tot(0"'') — C/fin(T'fin), where C/fin(T'fin) = C^tot(?fin, a) is the 
energy of the final equilibrium state. Because the entropy is conserved during variations of parameters, the optimal 
final equilibrium state will have a density matrix 

exp[-/5fi„Htot] f„ 
Pfin = , (7.7) 

involving the temperature Tfin — l//3fin, which is determined by constancy of the von Neuman entropy of the total 
system: 

StotiT) = Sb{T, 7 = 0) + SpiT, ao) = Stot{Tnn) = ^B(Tfi„, 7 = 0) -I- SpiTf^^, a) (7.8) 

with Sb{T,j = 0) = 7r^r/3?7,A from ( 3.18 ). Remembering that the level splitting A of the bath modes is very small, 
one can solve 

Tti,, = T+^[Sp{T,ao)~Sp{T,a)] (7.9) 

This yields 

|>V,„ax| = Wo + f/tot(T, ao) - f7tot(T, a) - TSp{T, ao) + Tp{T, a) 

= Wo + Ftot{T, ao) - Ftot{T, a) - Wo + Fp{T, ao) - Fp{T, a) (7.10) 

In the last step we canceled the contributions of the unperturbed bath. Not unexpectedly, the result depends on the 
free energies of the total system. 

Notice that for a cycle consisting of the changes ao — > a and (much) later a ^ ao, the maximally extractable 
work becomes the sum of the amounts of work Wo(ao — s- a) + Wo (a ao), so for cycles in principle all work can be 
recovered. 
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a. Values at high and low T 



In the classical limit the free energy is given in ( p.52| ), just the value for a harmonic oscillator, whether or not 
coupled to other harmonic oscillators. The maximally extractable work is 



\[-^ + Ml o.o)]T = -alT + 0{al) 
2 1 — ao 4 



(7.11) 



At low T the difference in free energy of the total system between the equilibrium states at the initial and final 
value of the spring constant is 



^tot(a) - ^tot(ao) = Fp{a) - Fp{ao) 



The energy added at t = is 

Wo = X - - 1 ao - — 
2 ao Zttj 



ha 
27r7 



ao , 1 + w 
— in 

w 1 — w 



2w 



■In 



1 + w , 
I — w' 



— In- 

w 1 — w 



4(1 



In 



1 + 



+ (ao + 2aa) 



2w 1 — w' 

So by making the instantaneous change in a the maximally extractable work ( [7.10 ) reads 



6ha 



ha 1 , 

-A J 1 

47r7 



1 — 1 + w , 

■ m ■ 



2w 



\ ~ w' 



6ha 



(7.12) 



(7.13) 



(7.14) 



If we let the system relax, this will run away in the bath on a timescale t^. By making clever subtraction schemes, 
we may recover some of it, and in principle all of it. 

For the case T = Tfin = Eq. ( |7l0| ) merely says that that all energy exceeding the groundstate energy of the new 
system can, in principle, be extracted. 



3. Work extraction by further sudden changes 



Here we present the formalism of work-extraction via sudden changes of a parameter. Besides presenting the 
general setup, we will display the validity of the Thomson's equilibrium formulation of the second law within the 
present situation. 

Let there be a closed system with a Hamiltonian 7i in a state p{ti) at the moment ti. Certain parameters of the 
Hamiltonian are varied in a very fast way such that for ti+ its Hamiltonian becomes Tii, but the state remains p{ti) 
due to the sudden character of the variation. The work done by an external source reads: 



Wi ^ tiip{h){Hi ~ H)] 



(7.15) 



In the second step the system is allowed to involve according to the new Hamiltonian Tii. At the moment t2 when 
the system reaches the state: 



p{t2) = e-*(*^-*i)"i/''p(ii)e*(*^-*i)^i/'' 
its parameters are suddenly returned to their original value. The work done in this step reads: 

The total work done by the source for this cyclic variation of the parameter reads: 



W = Wi + Wa = tr 



(7.16) 



(7.17) 



(7.18) 



where we have used Eq. ( 7.16| ). Notice that we consider the closed overall system, and only due to this fact the 
evolution of the system for times between ti and t2 is given by the Hamiltonian Tii. 

It is not difficult to see from Eq. (7.18) that the second law is satisfied for the present setup. Let us first assume 
that at the moment t the system was in the ground state of Ti.: p{ti) — |0)(0|. Then one has: 
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W = (0| e*fe-«i)-Hi A ji g-»(t2-ti)Hi/fi _ Ji |0) > 0, (7.19) 

just by the definition of tlie ground state. The same statement, namely W > 0, holds when p{t) is the Gibbs 
distribution of the initial state at positive temperature T = 1/(3: p{ti) = exp(— /37i)/Z, Z = trexp(— /37i) p3|| . 

Our work extraction mechanism involves a second change of the spring constant, which is cyclic: at time ^2 we 
impose a jump a — > 02 = a(l — 02) and it keeps that value, until at ^3 it is put back to a. The work involved in this 
cyclic two-step process is 

AW = Utotitt) - Utotit^) + (7tot(4) - Utotit^) = i(a2 - a)[{x^)t, ~ {x')t,\ (7.20) 

The change in particle energy between and tt, is 

AC/ = ^mt. - + ia[(x2),3 - ] (7.21) 

Thus the change in heat during the work extraction process is 

AQ = ^{^V^. - (P'>tj + \a2{(.x\., - i,x\,\ (7.22) 

The values of {x^) and (j?) at time t2 are set by the spring constants oq and a solely, and can be taken from 
previous section. When we take is large, we can take for that situation the limiting values for a system with spring 
constant 02- We then find 



Ai7 = irjaa) - -r„(a) + —TJa^) - -r^(a) - CeC^) (7.24) 

2 PV 2^ 2 ^ 2a2 ^ ' 2 "'^ ^ 27r7 ^^2m' ^ ' 



AQ = ^TM - -TJa) + -Ua,) - ^r.(a) - ^ CJ^:^) - ^ C.(^) (7.25) 



a. Classical regime 
In the classical case with strong damping we use ( |6.13 ) and obtain 

AW = ^alT - ^aoa2Te-^''''-/'^ (7.26) 

This has a minimum 

a; = AW,ni„ = -^a^T e-^"*^/-^ (7.27) 
z o 

Work can only be extracted under proper conditions, that is to say, when a2 is between and ao exp(— 20^2^2) ■ Else 
the cycle a ^ 02 ^ a disperses energy. 

The heat absorbed by the subsystem is at linear order in ao insensitive to the work extraction, 

AQ = AU = -^aoTe-'^"''/'' (7.28) 

There are two cases: 

• ao > 0: The central spring is stiffened, a > ao, and energy is supplied at i = 0. This energy leaks away, mostly 
as heat into the bath, (AQ < 0) and partly as work extracted from the total system {AW < 0), more precisely, from 
the particle. 

The ratio of extracted work to maximally extractable energy is, in the regime where ao is small but finite and 

T » Uluo, 

^^-^ = le--*^ (7.29) 



|W„,ax| 2 
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So our mechanism extracts maximally 50% of maximum; for doing this it must start immediately (^2 = 0) and last 
as long as possible (^3 = 00). 

• Q!o < 0: the central spring is weakened a,t t — 0. Energy is taken out from the system. The amoimt, ir(|ao| ~ciq) 
is less than the amount that could have been extracted by an adiabatic change, ^T(|ao| + ^Q^o)- After that has been 
done, heat flows from the bath to the particle (AQ > 0), as if the particle were at a lower temperature. In the course 
of this process work can be extracted, maximally the absolute value of Eq. ( [7.27 ), as is usual for two-temperature 
systems. The basic issue to extract work is to have a mechanism that, given the initial change in a, is the closer to 
adiabaticity. 

These conclusions also hold for moderate and weak damping. 



b. Low temperature regime 



At T = we get from Eq. ( |7.23D 



AW 



ha 
27r7 



Zm 



where Cx{0) is defined in ( |6.25 ) and Cx{t) in (6.6). The minimum occurs for 



"0 ^ (1*2. 



The change in heat and internal energy is 



AQ = AC/ 



7*2. 
27r7 2m 



(7.30) 



(7.31) 



(7.32) 



We consider again the separate cases: 

• ao > 0. The spring is stiffened at * = and energy Wo is supplied. At t = 0+ the energy is, to linear order in 
ao, equal to its final value, since Ce{0) — 0. This changes since energy comes from the bath, on a timescale Tp which 
is short for strong damping, where Tp = m/j. In the initial time regime 7*2/2771 < tq this mainly leaks away to the 
bath, and a small part can be extracted as work. In the regime tq < 7*2/2777 < cTmin this also happens, but the energy 
of the subsystem goes below its final value, so the particle becomes "too cold". In the final regime 7*2/2™ > fJmin 
energy flows back to the particle and again a small part can be extracted as work. This is then work extracted from 
the non-equilibrium bath, and the surprise is that this can be done although initially energy was put on the particle. 
This recovery of energy stored in the bath is a quantum effect. 



AW < means that work is extracted from the total system. Eq. ( |7.3lD is the maximally extractable amount of 
work with the pr esent mechanism. As efficiency factor way may normalize with respect to the maximally extractable 
energy from eq. ( 7.11 ), the energy that would otherwise leak away into the bath. 



-AWn 



(7.33) 



When one starts the extraction quickly after the initial change (ri « 0) one can still get half of the work back in this 
way, the same rate as in the classical regime. Even more is obtained when one still starts at *2 = but stops at the 
moment that the energy current goes no longer towards the particle, but away from it, i.e. at 012*3 = Umin- One then 
has 



AW = -^a2a[(a;^)t=0;ao - {x'^)t3;a2] = ^—Wl ~ aoa2[l - Cx{L^2t3)]} > -^^"ofl ^ C'2;(t^2<3)]^ 



The efficiency is 



-AW 



:[1-C,(u;2*3)]' 



(7.34) 



(7.35) 



which indeed has the proper behavior for *3 — > and *3 — > 00. The maximum is, due to (6.35|) and (3.3S), 
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?7max = ^[1 - C'4a,„i„)]2 = 0.5961 (7.36) 

the maximum exceeds the classical efficiency 77 ~ 1/2. So the quantum statistical excess energy flow from the particle 
into the bath indeed allows a more optimal recovery of energy initially put on the particle. 

The most interesting feature is that with the present mechanism it is also possible to extract work solely from the 
bath, a mechanism forbidden by the original Thomson formulation of the second law. Indeed, after 102^ = Cmin energy 
will flow back from the bath to the particle. By starting the extraction mechanism at loi^i = fmin and exploiting all 



times after this for the work extraction, the maximal efficiency (7.33) is 

r?max = ^C^'l^Tmin) = 0.00422262576 (7.37) 

In contrast to the classical case, it goes to a finite limit when ao 0. This occurs because in the quantum case 
the energy ~ Oq, with respect to which the extracted work has been normalized, is one order of magnitude smaller 
than the initially supplied energy Wq ~ ao- 

• ao < 0: The spring is weakened and energy Wo is extracted. Some energy can be extracted. For 102^2 < fmin it 
comes from the bath, but in the regime ^2^2 > fmin the particle has an excess energy, which then supplies the work. 

4. Work extraction by smooth changes of the spring constant 

Let us now consider the case where the spring constant a{t) = [1 — a{t)]a, is slightly changed (|Q;(i)| <C 1) in a 
smooth manner, starting from the equilibrium state a(— 00) — a. In appendix A we derive for the rate of work added 
to the system 



dW _ dWrov dn 
"dT ~ ~dt ^ d^ 

where 



dWrcv _ 7 da(T) 



dt 2 m dr 

is the adiabatic ('reversible' or 'recoverable') rate of work and 

dn ha da{T) f°° ^ ,^ jt 



iT. + |^a(.)C.(0) 



(7.38) 



(7.39) 



Jo ^'^"'^^ 



dt 47rm dr ./n 2m 
is the rate of energy dispersion. 

a. Completed changes 

Integrating over the full change one has 



T)a(T) (7.40) 



Wrev = -^(a/ - cc{)T, - ^C,(0)[a^ - c^] (7.41) 



where = a(—oo) and o;/ = q;(cxd), and 

n=— / dna'in) dra a'(r2)C,( In - rsl ) (7.42) 

For a full process (covering the whole region where a' ^ 0) 11 is non-negative, since it is an integral over a non-negative 
function, 

n = ^ dC\A{0\^^CM) (7.43) 

27 J -00 
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where the Laplace transform Cx was defined in Eq. (6.51); it can be verified that ^Cx{iC) is positive for all real C,. 
Furthermore 



^(0 = / 



(7.44) 



The positive energy dispersion for a completed, non-adiabatic cyclic change of system parameters is the Thomson 
formulation of the second law. We see that a positive dispersion also holds for non-cyclic but completed changes, as 
is known to occur on general grounds. 



We can check previous case a{t) = aod{—t), a'{T) = —aoS{T), for which Eq. (7.13) is at T = equivalent to 

1 ha 



Wo = >V,cv + n = -aoTx + —alCi°^ 



(7.45) 



6. Incomplete changes 



Let US now consider the temporal build up of this result in the regime of strong damping. Let a have the form 

a{T) = anM^t) (7-46) 

where a„i is a small amplitude and ft, is a bounded function {\h\ < 1), with h{—oo) = 0, h'{oo) = 0. If h{oo) ^ 
it could be an error function; if h(oo) = a Gaussian. 1/17 is the typical duration time of the change. Using Eqs. 
( |6l8pl9|j6^ we get from ( |aT0| ) 



dn 

'dt 



47rm 
12n 



h'{nt) 



f2Tr-fT\ 



I a 



ha 



(7.47) 



where c'ljP was defined in ( 3.49| ); for weak damping it equals 3TT^/am/8j, while for large damping it becomes unity. 
The integrated effect is 



m 



hila'^ 



\ ha J 



nt 



2a 



[h'{nt)f 



-h'{m)h"{nt) 



at 



(7.48) 



For a completed change the second and third term vanish, leading to a positive energy dispersion, in concordance 
with ( |03| ). It is seen that then the standard behavior 11 ~ applies when T is large enough or when f2 is small 
enough. However, in the quantum regime where the duration l/fi is smaller than the quantum timescale = ^/r, 
the last term in Eq. (7.4S) dominates, with a new behavior 11 il^. 

Another new quantum effect is that at low T and for typical times the second term in (7.47) is larger in magnitude 
than the other ones. This too occurs since the integral of Cx{t) vanishes at T = 0, and leads to new possibilities, 
that we discuss now. 



c. Work extraction from a smooth cyclic change 

One definition for a perpetuum mobile of the second kind is that there is a machine that performs a cycle in which 
it receives heat from a bath and converts this fully into work done on the surrounding. Additional requirements can 
occur; we shall discuss them in section |9. 7] . Here we analyze whether such a full energy conversion can be realized 
in our setup. The aim is thus to have a cyclic change of a system parameter with the properties 

AC/ = 0, AQ = -AW > (7.49) 



Including in (A. 6) also the contribution of the momentum, we have to linear order in a 

Uit)^^T,{a) + -Tx{a) + —[a{^)CE{0)- dr a'(^ - r)Cs(r) ] (7.50) 
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Let us choose for ^(t) a curve in the shape of a double bell, and consider the system at some time ti after the first 
peak, and compare it to the situation at a later time t2, after the second peak. We require that the spring constant 
has the same values at these instants, a{Ti) = a{T2), where ti^2 — 7^i.2/2m. This implies for the internal energy 

AU = U{h) - U{h) = ^ r dT[a'(Ti - r) - a'(T2 - t)]Ce{t) (7.51) 

For slow changes this can be expanded 

AC/ = ^{["'(^i) - ^'i^2)]cf + K(n) - a"(r2)]4'^ - i[a"'(ri) - a"'{T2)]cf} (7.52) 



where the coefficients are given in (6.53). 

Let us assume that we have two consecutive changes, characterized by a common bell-shaped function k{x) with 
|fc(a^)| < li but involving different rates of change ili.2, 

h{T) ^ k{nit) + k{Q.2{t - f^)) (7.53) 

where the parameters ili > > and ^2 > are such that the profiles have negligible overlap. Choosing the times 
ti,2 as 

we indeed satisfy the cyclic condition ot(T\) — oi(t2) = a^k(x). The difference in energy is 

AC/ = I f fc'(.) + I(^^i±^fc"(.) - l!M±i¥^^^'"(^)| (7.55) 

Yl-K \\ ha ) a I 

Let us assume that x is fixed but such that fc' < while k" > 0; this is possible because k has both convex and 
concave parts, implying that there is an interval with such behavior. In the example k(x) — exp(— a;^/2), there is an 
inflection point xn — 1 and one needs x > 1. Let us also assume that T <C ha/j. Then for 



Att^jT^ \k'{x)\ 

n\ k"{x) 



17i + f72 = (7.56) 



the first two terms can cancel, while the exposed correction term and the higher ones are small. Thus under these 
conditions it is possible to have a cyclic process with AC/ = to order Q;,n; if corrections are taken into account, 
the condition for cancellation is shifted an amount of order am and can again be met. Thus it is possible to start 
from the equilibrium state, make a first cyclic change of a and then a second, which process contains itself a cy clic 
change of a with At/ = 0. The work during this cycle comes solely from the energy dispersion. Using Eq. ( 7.48 ) we 
obtain the leading terms 



hoi 



12tt \ ha J { ' Js, ' " y_oo 2' " " k"{x) 



(7.57) 



The higher order terms are small for the same reason as above. If AW < this amounts to work exerted by the 
system on the environment. One can always have k"{x) small enough (by choosing x close to the inflection point) to 
make the combination of the Qi terms negative, and choose fl2 so small that the whole expression remains negative. 
So it is indeed possible to have a cycle where the extracted work |AyV| comes solely from the bath. 

We should stress that these work cycles are realizable only because in the first part of the process, for t < ti, energy 
was lost (yV(ii) — yV(— oo) > 0); a part of this is recovered. If, on the other hand, all the work is counted, then no 



work extraction is possible (W(oo) - W(-oo) > 0), c.f. Eq. ( 7.42 ). 

Alternatively one may conclude that there are non-equilibrium initial conditions (for instance, the state of the 
total system at time ti) which allow cycles which fully transform heat obtained from the bath into work done on the 
environment. 

From the analysis it is clear that under less strict conditions it is even possible to make a cycle that extracts work 
that comes partly from the bath and partly from the system ("efficiency larger than 100%"). 
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d. Perpetuum mobile with many work extraction cycles 



One can make several of these cycles. Even though the previous finding that complete cycles disperse energy should 
temper the hope to gain more work by doing more cycles, we consider the issue here, since work extraction from many 
cycles is one of the ways to express our unexpected results. 

Let there be Af cycles with a(r) = anih{T) having A/" non-overlapping bell-shaped parts, like in Eq. ( 7.53 ), where 
k{x) could be the Gaussian exp(— ix^). Each cycle is characterized its typical inverse duration time f2„ ^ a/7 and 
location f^, and each new cycles is slower than the previous one, r2„+i < fi„. For having cyclic behavior in the energy 
one finds from (7.55) that for low T and small but almost equal fin's, one should choose x close to the inflection point 
Xjf of k: 

-™^^^(^)^-^-<^) 

In order that this be small for all n one needs that fl/^ 3> jT'^/h'^a. Strictly speaking the a„ = amh{xn) are 
now not exactly equal; this can be healed by slightly adjusting the profile in the n'th cycle: k{nn(t — i„)) — > 
k{Q,n{t — t'^))k{xii) /k(xn), yielding a{Tn) = anik{xi[) for all n. For small Xn — Xif this correction factor is close to 
unity and can be omitted from the rest of the argument. 
Let us define 

T=^ (7.59) 
na 



Taking Eq. ( 7.48| ) at i = i„, where n cycles have been performed, and using that x — xa <^ 1 and that |f2„ — f2„+i| <C 
0„, brings for the yield of the n'th cycle 

n(t„) - n(t„_i) = \mnf^ + +1,-1^] (7.6O) 

127r 1^ 2a a^ J 

where 

/oo />oo 
d^WiOW h^&^\k'{x,,)]\ A^[k"{0? (7.61) 

00 «/ —00 

For having an equal yield per cycle, one demands 

n(i„) - n(t„_i) = Wcyeic = = ^^f\ (7.62) 

12777 12777 

where w > and w — T^v are dimensionless. It will turn out that there exists a consistent solution for v in some 
definite range. Assuming that Q{n) = fin is a smooth function of n one obtains 

2 

iinf^ + i2-nn' + i3^n^ = -w (7.63) 

a a^ 



-^-S/ , r^^. r .... (7-64) 



Solving for dn/dil and going to a new variable y = PMI/{2tt) one gets the total number of cycles 

h /-^^'"^/^^ Ayy 
f v + Iiy + hy^ 

The total yield is then 

mo.=AA>V.,.e = -^^f^/ d,— (7.65) 

127r 7 Ji3h^iu/27T v + hy + Isy-^ 

Here the minus sign indicates that work is performed by the system on the environment. This is possible because Eq. 
(7.48) expresses that, in order to make the work extraction cycles, one had to start from the equilibrium state and 
change a from «(— oo) = up to q;(ti). In this first part of the process energy was dispersed at an amount 
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(7.66) 



Notice that for v ^ {Phfti)'^ the extracted work becomes accordmg to (7.65) 

w,o,^-^^inl-nl) (7.67) 

so for rij^ <C ill there an almost perfect recovery, which is possible since the number of cycles is still large. 



For moderate v more cycles are possible, but less work is recovered. For the lower integration variable Eq. ( 7.58| ) 
gives 

For strong damping the physical timescale is = J /a. One assumes that fli is a large but finite number times I/t^. 
Choosing T <^ ha/j means that the upper integration limit yi = f3hCli/2TT is much larger than unity. But it is still 
possible to choose y^v <C 1, which is a useful condition for achieving many cycles. One then has for small v 

4\/7il3 iT 

which is indeed large. In the overdamped regime the yield can thus be expressed as 



mot — ———;== — VI — —J- I'-'Uj 

In the limit of weak damping we should notice that 

h = i2^fe (7.71) 

where e — amjy^ ^ 1 and where I2 is a numerical constant of order unity. Thus the work dispersed for achieving 
the non-equilibrium condition at t\ is 

n(t.) = +#^ (7.72) 

z47r cjQ 



where = ^Ja/m is the free oscillation frequency. Let us recall that = 2m/7 is the damping time. For T ^ hVli 
it holds that 

M=^=e^ (7.73) 

At the typical temperature T ^ hua this carries an additional large factor e. The yield per cycle carries a factor 1/e, 
so this total yield is independent of e. For small v it reads 

which can be comparable to the dispersed work, but it is always less. 

Summarizing this section, we have investigated the presence of many work extraction cycles both in the strongly 
and weakly damped regime. At low T their maximal number can be large but it is finite. The divergence M ~ 1/T is 
probably cut of at low enough T when the amplitude of our changes is small but finite. When more than M cycles 
are made, the possibility of work extraction is lost, because of the dispersion inherent to cycles. At moderate T the 
possibility of work extraction by cyclic changes is quickly lost; it is a strictly quantum effect. 



8. ON EXPERIMENTS TO TEST THE BREAKDOWN OF THE SECOND LAW 

In this section we will briefly comment on practical realizations of the low-temperature, non-weakly damped quan- 
tum Brownian motion. We do not intend to make detailed proposals for experimental setups, but we will mention 
certain fields, which according to commonly shared experimental views, display the above-mentioned strong-coupling 
and/or low-temperature regime. 
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1. Once more: the characteristic timescales 



Let us first recall once more that there are several important time-scales in the problem, tq is the characteristic 
time brought about by the external potential, which the particle will have if there is no interaction with thermal bath. 
For reasonably simple confining potentials there is only one such a time. In particular, for the harmonic external 
potential V{x) — ■^ax^ it is read tq = l/oJo = ^/m/a. Since no indications of damping are seen in this time, it can 
have a physical meaning only for very weak damping: 7 — > 0. 

If damping is large, then the characteristic dynamical times are Tp = m/j, Tx = 7/a- The overdamped regime 
appears with Tp <^ t^, and in this case Tp and t^ can be interpreted as the relaxation times of the momentum and 
coordinate, respectively. 

In contrast, very weak damping means 7 — *■ 0, and the damping time t^; ~ Tp = m/j is the longest characteristic 
time. For intermediate values of 7 the characteristic dynamical times in the overdamped regime are l/wi,2 defined in 
Eq. (3.31), and for the underdamped regime they are given in (3.38). 

The issue of this work is to consider the regime where another timescale, the characteristic quantum timescale 
Tn ~ fi/T, plays a dominant role. In particular, this timescale governs quantum correlations of the bath | |50| . The 
high-temperature classical case naturally corresponds to <C To,Tp,Tx and there quantum correlation effects can be 
neglected. 

Equilibrium quantum thermodynamics is recovered in the limit Tfi,TQ <C Tp, which means that the momentary 
motion of the Brownian particle practically does not notice damping, though it does so at long times. Obviously, 
this condition cannot be satisfied at low temperatures. In the present paper we are interested in the regime where 
both damping and quantum correlation effects are important: Tp ~ r^, ~ tjj, including possibilities of Tp <C Tx or 
Tp,Tx ^ Tfi, where our results are only strengthened. As noticed at the end of Section 7, new possibilities for work 
extraction has been found to occur in this last domain, where the inequality Tp <^ Tx Tn says that the noise is 
(anti-)correlated throughout the systems relaxation. It thus looks more like a quenched random variable than an 
annealed one, thus not at all behaving like a white noise, the standard ingredient needed to derive from a Langevin 
equation a Gibbsian equilibrium state. 

Before proceeding with concrete examples, let us just notice that there is nothing exotic in the quantum time-scale 
itself: Tn = 7.6 ps at T = 1 K, which is fully in range of the modern technologies. 



2. Possible experimental realizations 

a. Josephson junctions 



The first example to be discussed are Josephson junctions [ [L7yiq , |53| , |54[ |. This well-known phenomenon represents a 
standard example of quantum Brownian motion. The Josephson junction consists of two superconductors separated 
by a thin insulating barrier. Cooper pairs of electrons (or holes) are able to tunnel through this barrier, thereby 
maintaining phase coherence in the process, and leading to a possibility to have superconducting current. There is a 



direct map between properties of this junction, and the standard model of the quantum Brownian motion 1 17 lq , |53| , |54| . 
In particular, the coordinate x can be corresponded to the phase difference of the Cooper pair wave functions, the 
friction founds its place as resistance, mass is related to capacitance, and the current noise has the standard spectrum 



( ^.25 5.26| ), and can be related to ri{t). Under certain well-defined conditions one can neglect tunneling of the phase 



from one metastable state to another p^ , and consider it in a confining, nearly-harmonic potential. This system 
couples to the environment, which acts as the bath of our theory. In practice one can notice the occurrence of strong 
coupling at low T since then a careful shielding of the sample is needed in order to prevent an influence of the 
environment to the measuring apparatus. It appears that the non-weakly damped and low-temperature limits are 
well-known for Josephson junctions, and were a subject of rather long experimental activity |17 1^ , |5^ . For example. 



the following regime was explicitly realized as a condition of "really-quantum effects" |18| : Tp ~ 0.1 ps which is smaller 
than Tn at 1 K. The ratio Tp/Tp need not be of order one, but can vary significantly (from 0.1 to 10) depending on 
the construction of the junction; for details see [Q. In experiments reported in ||l^ the authors achieved Tn/Tx ~ 10 
at r = 1 K, and Tp/Tx ^ 0.1, which is a typical overdamped, low-temperature case. Notice that these experiments 
were among the first ones, where the spectrum of the low-temperature quantum noise was measured and found in the 
perfect agreement with the assumed standard form of the quantum Langevin equations. 
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3. Low temperature electrical circuits 



Experiments on mesoscopic, low-temperature electrical circuits |55 5^ provide yet another example, where non- 
gibbsian values of were clearly observed, and found in a good agreement with the theoretical predictions. We 
recall that the linear RLC can be mapped to the harmonic brownian particle: the coordinate x and the momentum 
p of the particle correspond to the charge and the current of the circuit, and m and a are directly connected with 
the inductance L and the inverse capacitance 1/C of the circuit (see also our discussion after Eq. (3.1)). Finally, the 
damping constant 7 corresponds to the ohmic resistance R. One notices that the (quasi) Ohmic limit, where T is the 
largest characteristic frequency of the problem, is conveniently realized in the present context. 

First of all we notice that for experiments described in Refs. |5^,Q all the relevant characteristic time-scales have 
basically similar values: ~ Tp ~ r^; = 10~* s, which makes the situation especially relevant for our purposes. 

Here we will briefly discuss the possibilities of experimental detection the Clausius inequality violation at low 
temperatures, since this seems to be the simplest possible issue. Moreover, the mo st ev ident situation is realized upon 
a slow variation of the inductivity (mass) L, where for T ^ — according to Eq. ( 4.25 ) and in the clear contrast with 
the Clausius inequality dQ < — one gets a finite positive heat provided that dm = dL > 0. One needs to observe 



^) and {p^) for several different values of the inductivity (mass) L. This is s ufficie nt to recover the corresponding 
changes of the average energy, as well as to recover the work according to Eq. ( 4.24 ). The heat is then obtained by 
subtracting the work from the energy. In the second step one can check the consistency of the results by observing 
directly the work done by the external source. Altogether, the challenge of the main experimental observation is in 
observation of the variances. 

In Ref. psf the authors considered mesoscopic electrical circuits in the context of single charge tunneling. The 
used circuits had thickness of the order 10 nm and wideness of the order 1 /im. The observations allowed indirect 
determination of (x^). With the subsequent improvement made in the correspondence with the theoretical 
expression (3.47) is perfect. The observations were done with C — 1/a =4.5 fF, L =4.5 nH and for i? = 7 in the range 
10^ — 10"^ kfi. For damped circuits the relative importance of damping is quantified by the quality factor Tp/ro, which 
in the above range of parameters varies from 10~^ to 10~^. To avoid thermal noises the circuits were cooled up to 20 
mK. At such a low temperature quantum effects are dominating, since the quantum time-scale th = h/T ^ 10^^ s is 
larger than the other ones, tq ~ 10~^ — 10"^*^, Tp ~ 10^® and ~ 10~^ s. To get an idea for the magnitude of the 
expected effect, let us estimate the outcome for AQ ~ L dQ/dL. With the above parameters and i? = 7 = 10'^ kfi 
one gets from ( |4.25| ) AQ - 10" J= 1 eV. 

Since for mesoscopic circuits the formula for was already verified, it is now a matter to perform three measure- 
ments (the equivalents of Tp, and the work production) on a single sample, to verify unambiguously the breakdown 
of the Clausius inequality. 



a. Trapped ions 

As another, more elementary example one can mention a trapped ion immersed in a photon bath. Taking as an 



estimate the mass of proton (to = 10 kg), and 7 = 10 kg s ^ |12 , one gets ^ t^j at 1 K, so the quantum 



coherence effects are still active. The ideal example of a harmonic Brownian particle will be an ion trapped in a 
so-called Paul trap Q| , or an electron or ion in a Penning trap |4|] . These electromagnetic traps are nowadays well 
realizable and suited for variation of parameters. In particular, high quantum number Rydberg states may have a 
long lifetime and a strong coupling to the radiation field. 



9. ON THE FOUNDATIONS OF THERMODYNAMICS AND PERPETUUM MOBILE 

This section intends to summarize to what extent the standard relations and laws of thermodynamics can be applied 
to a quantum Brownian particle. There are many formulations of the second law, and some of them have been found 
invalid in previous discussion. One may go to the extreme limit as saying that there is no motivation to discuss a 
thermodynamics in the way we did. To show that there is justice in doing it, we summarize our results in the light 
of standard thermodynamic wisdom, and point at the agreements and contradictions. 

For a general, pedagogic text on the history and today's status of thermodynamics and the second law, we refer to 



the recent work by Uffink |Q. For a collection and discussion of the original papers, see the book by Kestin |47|. A 
very recent discussion on the base of the axiomatic thermodynamics was presented by Lieb and Yngvason |58[ . For a 
discussion of what can be meant by "the" entropy of a system, see [p9i . 
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1. Has thermodynamics been violated or did it never apply? 

The conclusion of our analysis is that thermodynamics does not work when, in the quantum regime, ones considers 
the Brownian particle in its reduced Hilbert space, thus summing out the bath variables of the total system. This 
implies that the characteristics of the particle are directly observable which is indeed the case with the standard 
examples of the brownian motion. There are, however, situations, where only some composite (system-plus-bath) 
quantities are measured, and the need for a separation between particle and bath is questionable (there still can be a 
possibility that such a separation can be done on a different, more coarse-grained description of the overall system, but 
we will not enter into that discussion here). This is the case with a Kondo system, where the measured quantity is for 
instance the magnetization, which is set by the magnetic impurity and the bath together (i.e. it lives in the common 
space of the particle and the bath). Also for the dressing of a 'bare' electron by photons, it is the standard practice of 
the quantum field theory to consider the dressed mass and charge as directly observables quantities. However, when 
the system is a Josephson junction or a mesoscopic circuit, its own characteristics are perfectly measurable, so there 
is an important case to make. When looking at the budget of the junction alone, one has to keep in mind that it can 
exchange energy with its environment. At low enough temperature this mechanism displays unexpected behavior and 
is responsible for non-thermodynamic characteristics. 

If one accepts to consider the subsystem as a Brownian particle immersed in a heat bath, then first it should be 
noticed that the particle acquires a cloud of bath modes around it. This dressing is a manifestation of the (strong) 
damping of the particle by the bath. One can then ask the question: "if thermodynamics does not apply, where was 
it lost?" If no technical errors have been made in our derivations, then the answer must be: "It never applied" . 

An argument in favor of this point of view is the fact that the interaction energy is non- vanishing, thus violating a 
well known condition for the derivation of the standard equilibrium thermodynamics. The interaction energy reads 

c/int - t/p - - ir, (9.1) 



The result is given explicitly by Eqs. (3.45) and ( 3.48| ). In a strict formulation of thermodynamics one requires that 



the equilibrium value of the interaction energy Uint is negligible, in order to separate what is meant by the system 
from what is meant by the bath. In our case this would imply 7 — s- or T — > oo, and indeed in both limits standard 
thermodynamics is recovered. However, in general the same system has a non-zero interaction energy and, in the 
strict formulation, should not be considered to be thermodynamic at finite 7 or T. Whereas the limit of large T 
can be naturally achieved in practice, the weak-coupling limit 7 — > is much more difficult to realize, since coupling 
constants are generally fixed numbers whose magnitude cannot be manipulated at will. Likewise, all physical systems 
having a non- vanishing interaction energy with their baths should then not be considered as thermodynamic systems. 
This would apply to very many systems at low enough temperature, leaving an uncomfortable situation with respect 
to the well behaved high temperature properties of the same systems. Already at high T the analysis of the Clausius 



inequality in section 4. 2 c gave a compelling argument in favor of our choice for the Hamiltonian of the subsystem, 
leaving the rest of the total Hamiltonian for the bath. 

Let us inspect in some detail the weak damping limit 7 <^ ^Jam. Here it can be shown that, even at T = 0, 



one has C/int ^ U ^ since, due to (3.48), C/int ~ 7 but U implying that the condition for the application of 

thermodynamics is almost fulfilled. Nevertheless, the Clausius inequality is typically violated, by an amount of, again, 
order 7. This argument somewhat weakens the point of view that thermodynamics should never apply. One can also 
argue that it is compatible with it, since the violation of thermodynamics is of order of the small interaction energy. 

Putting all arguments together, we reach the unavoidable conclusion: there are principle problems to define ther- 
modynamics at not very large temperatures and, in particular, in the regime of quantum entanglement. There is no 
resolution to this, and thermodynamically unexpected energy flows appear to be possible. 



2. Zeroth law 



The zeroth law is often said to state that in an equilibrium situation there will be a unique temperature. A 
standard formulation is that if two bodies are each in equilibrium with a third body, then they are also in equilibrium 
with each other, and the three bodies have a common temperature. Let us look, however, at a careful formulation 
out of equilibrium: // two parts of the system have an infinitesimally small temperature difference, then they will 
spontaneously equilibrate and reach a common temperature. 

For the (nearly) harmonic situation two different effective temperatures Tp and can be related to the momentum 
and the coordinate. Recall that these temperatures arise from the generalized form (4. 14 -4. IS) of the Clausius 
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inequality. The legitimation of such a definition of effective temperatures is confirmed by their successful use in glassy 
thermodynamics ]23|,p5t. 

In our case, Eg. ( |3.5^ ) shows that Tp deviates at large T from by a term with A — fi'^ jVlm. So for any 

infinitesimal e, the regime T > AT jt indeed has temperatures Tp and Tx that differ less than e. However, since they 
are parameters of the steady state, they will not equalize spontaneously, in conflict with the above formulation of the 
zeroth law. Instead, they become more and more different from each other at lower temperatures, and at zero bath 
temperature Tp and are both finite but different from each other. The fact that they remain finite just indicates 
that the corresponding quantum state does not have sharp values for p and x\ so this is a consequence of quantum 
complementarity. The fact that these effective temperatures take non-gibbsian values is a consequence of the quantum 
entanglement. In the Gibbsian limit of weak coupling (i.e. 7 — *■ 0) for the harmonic oscillator, both temperatures T^^ 
Tp tend to their common value ^?iwocotanh(ifi./3wo) of the harmonic oscillator weakly coupled to its bath. 

We should mention that the existence of the zero law is frequently viewed just as an axiom, but under certain 
conditions it can be derived from the second law (the entropy of a closed system never decreases) § . As we mentioned 
already, this derivation is based on the use of a weak interaction between the particle and its thermal bath. It confirms 
that if this weak-coupling condition is valid, then the two effective temperatures are indeed approximately equal. 

The fact that these unequal effective temperatures do not cause heat currents that equalize them, as would be 
required by standard thermodynamics, is reminiscent of the classical paradox that atoms should radiate, but, being 
in the quantum regime, they do not. 



3. First law 



The first law relates the change of system's energy into the heat added to it and the work done on it. It can not be 
broken, since it is a direct consequence of energy conservation, a central concept in quantum mechanics. Nevertheless 
the formulation of this law is not merely a tautology, because it allows to separate clearly those ingredients of the 
energy change, which arise from non-observable degrees of freedom (heat obtained by the Brownian particle from the 
thermal bath) and external sources (work done by them on the whole system) . Our identification of the energy of the 
subsystem as the expectation value of the Hamiltonian Ti. was supported in section 4. 2 e by requiring application of 
standard thermodynamics at high T, and is imposed by the form of the Wigner function. We stress that, given this 
identification of energy, our identification of the heat dQ added to the subsystem and the work dW done on it, are 
well accepted and widely discussed in literature, see e.g. the books by Keizer p%, Balian |M and Klimontovich M. 



4. Second law 



Let us stress that there are many formulations of the second law. There are several formulations of the second law 
that are, at least apparently, violated by the solution of our problem. 



a. Thomson's formulation, Kelvin's principle 

The formulation by William Thomson, the later Lord Kelvin of Largs, is; It is impossible to perform a cyclic process 
with no other result than that heat is absorbed from a reservoir, and work is performed. An earlier and more particular 
version of this statement appeared due to Carnot. 

For general quantum systems starting from the equilibrium state, this can be proven mathematically ; a simplified 
proof will be presented elsewhere ||60(|. In our setups it can always be verified, see for example the fact that the energy 



dispersion ( 7.43 ) is non- negative. After finishing the cycle the bath is not exactly in its Gibbsian state, but it is still 
very close to it, because the bath is extensive. Basically the dispersed energy has run away to infinity, leaving the 
system locally again in a Gibbsian state. This implies that also successive cycles will always disperse energy. 

However, out of equilibrium Thomson's formulation appears to be endangered. The first point to notice is that this 
can already occur at the classical level. The reason is simple. Consider, as we did in section [7. 3| , a sudden weakening 
of the central spring. In doing so, energy is extracted from the system, but, due to the sudden nature, it is not the 
optimal amount. One can improve on this by making the following cycle: quickly put the spring back at its original 
value, and then make the change slower. This cycle, that started in a non-equilibrium state, will yield work, and this 
work comes from the bath. We conclude that the Thomson formulation can only refer to system changes on long 
enough timescales, such that the initial state is practically in equilibrium. 
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A more exciting violation of the non-equilibrium Thomson formulation was observed for smooth changes of the 
spring constant at low enough temperature. In section 7.4d we discussed the case of A/" 3> 1 bell-shaped cycles in 
the spring constant; each cycle has two inflection points, pre-peak and post-peak. Starting in the Gibbsian state, the 
first cycle up to the post-peak inflection point is considered as a mechanism that produces for us a proper 'initial' 
non-equilibrium state. If the typical duration of the successive cycles increases, parameters can be chosen such that 
after each return of the spring constant to its post-peak inflection point value, the system has the same energy, while a 
prescribed, flxed amount of work is extracted. There can be A/" ~ 1/T of these cycles, which can be large at low enough 
T . They extract heat from the bath and convert it fully into work, forbidden by the general (i.e. non-equilibrium) 
Thomson formulation. It could be checked that the total amount of extracted work is less than the energy dispersed 
in the first part of the first cycle, so energy conservation is not endangered. The interesting fact is nevertheless that 
there can he M ^ 1/T of these cycles, which can be large at low enough T. Actually, making more cycles implies a 
smaller total extracted work ^ l/A/"^, since these cycles themselves lead to additional dispersion. 

In contrast to the violation of the Clausius inequality, to be discussed below, the violation of Thomson's formulation 
is a consequence of both quantum regime (low temperatures) and the non-equilibrium character of the whole system 
(particle and bath). Indeed, any amount of work extracted by means of the particle is in fact extracted from the 
whole system. If this global system is in equilibrium (namely it is exactly described by the Gibbs distribution) , there 
will be no possibility to extract work by making a cyclic change of a system parameter, since this formulation of the 
second law applies as well to any closed equilibrium system [^,|o) . On the other hand, the full account of quantum 
effects is necessary to show our work extraction, since it disappears in the Gibbsian limit, namely both at moderate 
temperatures and/or for weak coupling to the thermal bath. 



b. Clausius ' Principle 

Clausius states: It is impossible to perform a cyclic process which has no other result than that heat is absorbed 
from a reservoir with a low temperature and emitted into a reservoir with a higher temperature. 

The gained work could be used to drive some frictional process at a higher temperature, which would turn it into 
heat, in conflict with this principle. Nevertheless, this principle is obeyed at high temperatures, and only violated in 
the quantum regime at low T. 



c. Clausius inequality 

This formulation claims that in any thermodynamical process (in particular, for variation of a system parameter) 
the amount of heat received from the thermal bath by the particle is limited from above by the bath temperature 
times the change of the von Neumann entropy of the particle: 

dg < Tds.,N- 

A particular formulation of this law is that no heat can be extracted from a zero-temperature thermal bath, it can 
only be dumped in it (i.e., then it is impossible to have dQ ^ 0). This situation is particularly interesting, since it 
does not employ in any way the concept of entropy, and therefore can be applied to situations, where entropy is not 
known, or not well-defined. Physically it is also easy to understand. The energy of the cloud of bath modes around 
the subsystem will change if a system parameter is changed, even at T = 0. This change in energy of non-observable 
modes is heat, and it can be positive or negative, depending on the sign of the change. In one of the cases energy 
from the cloud will increase the subsystem's energy, violating the Clausius inequality. For a closed system dQ goes 
to zero, and one recovers from the above inequality the more standard formulation dSvN > 0, which appears to be a 
particular case of the Clausius inequality. 

We have shown that both those general (all T) and particular {T = 0) formulations are violated in quantum 
case. Although at high temperatures these violations are small, they nevertheless do exist. For researchers who are 
reluctant to follow our identification of the effective temperatures, it will perhaps be hard to agree on the violation 
of the zeroth law, discussed above. However, the violation of the second law, which also sets in at arbitrarily 
large temperatures, should be easier to accept, since the Clausius inequality does not employ the notion of effective 



temperatures. Moreover, in section 4. 2 d we have discussed a formulation which compares only equilibrium systems. 

For our harmonic system we succeeded in generalizing the Clausius inequality, involving two temperatures and two 
entropies, in the very same way it was done for glassy systems and which applies to black holes. 
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In hindsight, the derivation of the Clausius inequahty is nontrivial in the case under consideration. In standard 
thermodynamics one formulation of the second law is that the total entropy of a closed system cannot decrease. 
When applied to a subsystem coupled to its equilibrium bath, this immediately leads to the Clausius inequality. 
Hereto one makes two assumptions: equilibration of the bath and additivity of the entropy. Let us follow the 
subsequent steps. Because of its equilibrium nature, the heat received by the bath is associated with an entropy 
change, dQbath — TdSbath- Energy conservation says that S^Qbath +dQ — 0, where the latter is the heat received by 
the subsystem. This implies 

dS-^^dS + dSbath = d^tot > (9.2) 

In the world of quantum entanglement, however, both assumptions are less obvious. First, it does not hold that 
Stot = S ^ Sb- We have shown this explicitly, since at T = one has S'tot = S's = 0, but S — S^n > 0. Both the 
fact that energy is not quickly redistributed in the bath and the non-additivity of the entropy imply that there is no 
a priori reason to expect that the Clausius inequality is satisfied. In concordance with that, we have shown that it is 
indeed not valid. 



d. The rate of energy dispersion is non-negative 



In section 7. 4 b we have pointed out that at low enough temperatures the rate of energy dispersion can easily be 
negative. This holds even when one starts in equilibrium. Thus non-negativity of energy dispersion cannot serve as 
a universal formulation of the second law. 



e. The total entropy of a closed system cannot decrease 

The most standard formulation of the second law is that the (coarse grained) entropy of a closed system cannot 
decrease. In classical physics for a subsystem in contact with a heat bath the equivalent is that the rate of entropy 
production of the subsystem is nonnegative. For situations close to equilibrium it can often be expressed as a bilinear 
expression in generalized currents, and the matrix elements are called Onsager coefficients; this matrix is positive 
definite in all known examples. 

The rate of production of Boltzmann entropy was also considered by us. In the case of weak damping there occur 
oscillations in the production rate around zero in each period; this sets in at moderate temperatures, and is akin 



to the oscillations in the energy, that occur already at any non-infinite temperature. In sections 3. 2d and 3.2e we 
have pointed out that even at low temperatures and in the limit of strong damping the rate of Boltzmann entropy 
production can be negative. So also this criterion does not qualify as a solid definition of the second law. 

We should stress that we did not find sensible production rates for other entropies. Perhaps not accidentally, the 
Boltzmann entropies of the coordinate and momentum sectors are the ones that enter our generalized, two temperature 
version of the Clausius inequality. 

In our setup the von Neumann entropy for the full closed system (fine grained entropy) should not be confused 
with the von Neumann or Boltzmann entropies of the subsystem, which pertain to the Brownian particle only. The 
von Neumann entropy of the full system is not altered by changing the strength of the spring constant. This entropy 
remains constant during the overall unitary evolution of the whole system, and also remains constant during variations 
of a parameter, since also there the overall evolution is still unitary. The formulation of the second law in terms of 
non-decrease of entropy definitely refers to the coarse grained entropy. In the classical situation the fine-grained 



entropy is conserved as well, by the Liouville dynamics. For more definitions of entropy, see 59 1. 

In passing we notice that if one starts from a Gibbsian state of the total system (central particle coupled to the bath), 
and changes a system parameter, then the conservation of entropy prevents the system to relax to a new Gibbsian 
state of the total system, since our total system is isolated. Nevertheless, the subsystem (the central particle) does 
relax to a state characterized by the parameters, which can be coded in the effective temperatures, of that would-be 
global Gibbsian state. It is the finite amount of energy dumped in the extensive bath that does not relax, since our 
bath lacks anharmonic interactions, or coupling to an external superbath. In contrast to a superbath, anharmonic 
interactions do not change the essence of the argument on the overall unitary evolution, conservation of both the von 
Neumann entropy and the energy. However, they can widen the set of observables having would-be Gibbsian values. 
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5. Third law 



This law claims that if the ground state of the Brownian particle is non-degenerate, then its von Neumann entropy 
is equal to zero. This is a direct consequence of the quantum Gibbs distribution, which predicts the pure vacuum 
state at low temperatures. In our case neither the von Neumann entropy nor the Boltzmann entropy vanishes when 
the bath temperature is zero. This occurs because of the quantum Brownian particle is in an entangled mixed state, 
and therefore cannot have vanishing von Neumann entropy. 

The third law is recovered when taking the weak-coupling limit. In that case Tp = — i7iti;ocotanh(i/3?ia;o)j 
implying that the parameter v of Eq. (4.34) takes the value v — i cotanh i/JTiwoj which causes the von Neumann 
entropy of the particle (4.35) to vanish at T = 0. In a certain sense the violation of the third law reported here for 
non-weak coupling is the most straightforward consequence of quantum entanglement. 



6. Perpetuum mobile of the first kind 



Taken literally, a perpetuum mobile performs perpetual, i.e., everlasting, motion. Nevertheless, rotational currents 
in ordinary superconductors, which may exist several days, are rarely connected to perpetual motion. We shall 
therefore employ the word "perpetuum mobile" for any principle that yields work. 

One speaks of a perpetuum mobile of the first kind when the first law is violated, leading to an everlasting 
performance of work without any cost. Such a setup is impossible in quantum mechanics, since it satisfies the first 
law by principle. So here is no issue in the question what "perpetuum" means precisely. 



7. Perpetuum mobile of the second kind 



Another formulation of the second law is: It is impossible to construct an engine which will work in a complete 
cycle, and convert all heat it absorbs from a reservoir, into mechanical work [ ]6]| . A machine which would do so is 
called a perpetuum mobile of the second kind, and the second law states that such a machine is impossible. 



a. "Perpetuum" mobile or perpetuum mobile? 



When the first law is respected but the second is violated, one speaks of a "perpetuum mobile of the second kind" . 
However, we wish to make some remarks on the word 'perpetuum'. Surely, in the eighteenth century such a perpetuum 
mobile was imagined, for instance, to cross the Atlantic by boat using only the energy stored in the ocean water. As 
such, there would be a basically infinite bath, and the mobile, if realized, could function perpetual, i.e. "for ever" . In 
general, when the bath is finite, it obviously has a finite energy at its disposal. In many setups, such as those with a 
finite rate of energy extraction, this implies also finite duration of the process. Thus even in the classical situation, 
the term "perpetuum" need not be a precise adjective for this type of mobile, and the point of view could be taken 
that a perpetuum mobile of the second kind need not function arbitrarily long, but must only work many cycles. In 
view of the failure to find so far any practical realization, this stretch of the definition seems allowable to us. 

In quantum physics the situation is even more clear. In a closed system the energy can never go below the ground 
state energy, so the amount of extracted work is always finite. As a result, there can never be an infinite amount 
of cycles for which a definite amount of work is extracted per cycle. One of the formulations of the second law is: 
"Perpetuum mobile of the second kind do not exist" . Such a strong physical statement must, of course, be richer than 
the general statement on the existence of a ground state. The crux is that already one of such cycles, that extracts 
work from a thermal bath, is forbidden. So, already in general, perpetual motion of the second kind does not have its 
literal meaning of everlasting motion; rather, it is a notion for a work extraction principle, and one cycle is therefore 
good enough for the birth of a 'baby perpetuum mobile'. 



b. The present situation 



For our aim the allowance of non-eternal duration of perpetuum mobile is relevant, since our effects only hold as 
long as both the particle does not relax, which happens on the timescale t^, and are quantum coherent, which involves 
the quantum timescale h/T. 
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We have discussed a work extraction mechanism that cychcly changes the spring constant in a certain time interval 
Each of these cycles is slower than the relaxation time of the system. When the quantum timescale Tn = h/T is also 
slower than the relaxation time, there occurs unexpected behavior: the contribution to the rate of dispersion inversely 
proportional to the duration of the cycle, normally the leading term, has a small prefactor quadratic in temperature. 
Therefore quadratic and cubic terms in the inverse duration also play a role. Out of equilibrium cycles have been 
designed where a constant amount of heat is extracted from the bath is fully converted into work, while the energy 
of the subsystem is at the end of each cycle back at the value of the beginning. In this sense, systems described by 
our models, with parameters in the appropriate regime, present at low temperatures true realizations of perpetuum 
mobile of the second kind. Probably, it is also possible to extract work both from the bath and from the subsystem 
("efficiency larger than 100%"). 

In a more stringent definition of perpetuum mobile one requires that in the cycles work is extracted "without any 
further change" . For our system this can be expressed as the requirement that the Wigner function of the subsystem 
be back at its original value. This would imply the requirements that AK = 0, AV = and AV = over each cycle, 
rather than only having AK + AV = 0, where K is the expectation of the kinetic energy and V of the potential 
energy. The question whether this extended constraint can satisfied by changing the spring constant and possibly 
also the mass, is left for the future. 

We should stress the conceptual difference between the present situation and the well-known case, where work can 
be extracted due to a temperature difference between two thermal baths [||. The latter is the standard setup for the 
thermodynamic heat engine: two baths are explicitly separated from each other, and therefore the whole system is in 
a nonequilibrium state, and can be used to perform work. If those baths are kept in a direct contact for a sufficiently 
long time, then they will go to equilibrium and after that no work can be extracted any more [^8| . In contrast, here 
we are have presented a case with one single thermal bath. After a sudden increase of the strength of the central 
spring, the central particle will go to equilibrium after one relaxation time or r^. But in doing so, more than 
the initial surplus energy is transfered to the bath, and in particular to the cloud of bath modes in its immediate 
surrounding. After a certain moment, this heat flow towards the bath stops, and then a smaller backflow occurs from 
the bath to the particle, before the whole comes into equilibrium. This backflow of heat is the mechanism that makes 
it possible to extract work from the bath by manipulating the particle, in a situation where this would be impossible 
classically. In particular, for smooth changes at zero temperatures the integral of this relaxation function is needed, 
but it appears to vanish, leading to a variety of new effects at low temperatures. 

On a thermodynamic level, the analogy with the classical case was strengthened because we could identify effective 
temperatures, though we also stressed that by themselves they do not tend to become equal to each other. All these 
intriguing aspects arise due to quantum effects, since we showed in detail that the same system coupled to the same 
bath displays at high temperatures the fully expected thermodynamical behavior. 

8. Perpetuum mobile of the third kind 

One can define a perpetuum mobile of the third kind when work is performed at the cost of a diminishing, but still 
non-vanishing zero-temperature entropy. This can, in principle occur in systems, such as glasses, which are able to 
relax to equilibrium, but are temporarily stuck in certain metastable states. Then the zero-temperature entropy can 
be used as a measure of this metastability [|9| . 

One could wonder whether our extraction of work is due to the present non-vanishing zero-point entropy. However, 
this is not the case, since for the purely Gibbsian case of particle and bath, the particle would have the same zero point 
entropy, but no work could be extracted. Moreover, in our case the zero point entropy is an indication of quantum 
entanglement and not of metastability. 

10. CONCLUSION 

This paper is devoted to the statistical thermodynamics of the quantum Brownian motion. The high-temperature 
case of this model can serve as a convenient pedagogic example, where almost all main statements of statistical 
thermodynamics are derived exclusively from the first principles. Among other advantages, such an approach makes 
possible to reveal the conceptual restrictions and limitations of the standard thermodynamical wisdom. With this 
aim in mind we focussed in the paper on the low-temperature (quantum) situation of the Brownian motion model. 

The stationary state of a quantum Brownian particle non- weakly interacting with its thermal bath is non-Gibbsian. 
It is this property which makes the quantum Brownian motion a challenging problem, and classical thermodynamical 



57 



wisdom appears to be inadequate even if the total state of the system and the bath is Gibbsian. Both the classical and 
the quantum Gibbsian thermodynamic theories emerge as particular limits in this more general setup. The classical 
Gibbs distribution with all its thermodynamic consequences is recovered for high temperatures, and the quantum 
Gibbs distribution is obtained for very weak damping. 

In section ^ we start from the quantum Langevin equation. At low temperatures this equation contains a colored 
Gaussian noise; because of quantum coherence, the bath cannot generate white noise even in the limit where the 
friction has no memory. To achieve this interesting situation, no more is needed than the observation that quantum 
mechanics applies (see detailed explanations after Eq. ( |).26[ )). Moreover, the quantum fluctuation-dissipation theorem 
predicts different time-scales of noise and dissipation at low enough temperatures ||J|,|l^,|ri| . This is how quantum 
coherence enters into the considerations. Its characteristic timescale is fi/T , where we had set the Boltzmann constant 
Ub equal to unity so far. Restoring it, we have 1% = h/{kBT) = (6.23/r) Kps. At T = IK one has tfi = 6.23 ps, in 
the range of typical microscopic processes in condensed matter. 

Since the stationary distribution is non-Gibbsian there are concept ual differ ences compared with the equilibrium 
case. For the harmonic potential one can define effective temperatures (4T, 4T) for momentum and coordinate. Both 
temperatures are different from the b ath tempe rature T. Generalized thermodynamic relations can be introduced, 
which take a quasi-Gibbsian form, Eq. ( 1.21 - 4.22 ), and are closely related to the ones in the thermodynamics of glasses 

MM; 

The inapplicability of standard thermodynamics is most clearly illustrated by the violation of the Clausius relation: 
6.Q < TdS: heat received by t he p article from the bath is restricted by temperature of the bath times the change of the 
particle's entropy. In section |4. 2| we construct an explicit example which at low temperatures realizes dQ > 0. This 
violation is significant at low temperatures, where quantum effects are relevant, and is small for high temperatures. 
It is important to notice that this violation exists already for the totally equilibrium (Gibbsian) state of the overall 
(particle plus bath) system (see our discussion after Eq. ( 4.2?! )). Since Thomson's formulation of the second law is 
valid for such a state |Q, we have the explicit counterexample showing that the very equivalence between different 
formulations of the second law is broken at low temperatures. 

For the dynamical consideration we start from a non-equilibrium initial state obtained from the total Gibbsian by 
changing the width of the confining potential. This change involves a small, controllable energy input, and can be more 
realistic than the hitherto studied case where particle and bath are initially uncorrelated. After the non-equilibrium 
state has been prepared, the ensuing relaxation of the particle presents a number of thermodynamical anomalies at 
low temperatures of the bath. First, energy put into the bath does not completely dissipates there (in contrast to the 
classical situation), and thus work-extraction from a single thermal bath is possible. This violates the second law in 
Thomson's formulation as applied to non-equilibrium. As a consequence of this, adiabatic changes of parameters are 
not the most optimal ones anymore. It is interesting to notice that at low temperatures a large (but finite) number 
of work-extracting cycles is possible. According to our opinion, this explicitly realizes the basic non-trivial content 
of the perpetuum mobile of the second kind, because any possibility for infinite number of such cycles is ruled out 
merely by the existence of the ground state for the overall system. 

The second aspect of the low-temperature relaxation is that no H-theorem exists at low temperatures, i.e. the 
properly defined entropy production appears to be negative for some times. This holds even in the moderately 
overdamped regime, when the brownian particle relaxes monotonously. (Recall that without any bath those frequencies 
are purely imaginary which leads to the known oscillatory behavior.) Within the underdamped (weakly-coupled) 
situation negative entropy production persists up to high temperatures, and disappears only in the explicitly classical 
limit. 

Let us recall that Thomson's formulation of the second law in its most general and universally applicable form 
In cyclic processes no work can be extracted from a closed equilibrium system, remains satisfied. This statement was 
derived in Ref . , and we will present a very simply proof elsewhere |]60| . It obviously applies to the analysis of this 
paper, since we start from a Gibbsian with modified spring constant, and the work extraction disappear when the 
change in the spring constant vanishes. Also for cyclic smooth changes that start from equilibrium we could verify 
the non-negativity of the energy dispersion. 

We now make some remarks concerning the definition of the thermal bath in our problem. The harmonic oscillator 
bath model, which was used by us, is technically convenient, but at the same time it possesses all relevant properties of 
a thermal bath, which are typically postulated in the statistical thermodynamics. The main of them is that the bath 
should have infinite amount of degrees of freedom, a necessary condition to ensure relaxation of the Brownian particle. 
On the other hand, the quantum Langevin equation, which is the starting point of our analysis, can be derived from 
rather different schemes (see e.g. pT|), since in a sense they are more universal than the detailed properties of the 
considered thermal bath. 

Finally, let us relax the conditions under which our results have been derived. We have already mentioned that they 
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hold as well for A'^ ^ 1 Brownian particles in an external potential. Though mutual interactions would complicate the 
analysis, it would not modify our basic statements. This can already be seen from the case of non-interacting harmon- 
ically bound Brownian particles: under a change of variables they become interacting ones, while the characteristics 
of the bath remains basically unaffected because it has very many degrees of freedom. 

Our findings on the non-thermodynamical character of the low-temperature brownian motion may have a wide scope 
of applications such as cooling, energy storage, and thermodynamical limits of low-temperature computing. Indeed, 
in the domain of information theory there is a large literature based on the fact that only the erasure of information 
must necessarily involve a dissipation of heat (see extensive reviews in 32 1 ) , and the claim that the erasure of one bit 
of information costs at least an amount of energy ksT \n2, the so-called Landauer bound. It is well known that this 
bound is based on a straightforward application of the Clausius inequality. From our observations it is clear that the 
Landauer bound can also be broken, and strong effects may occur at low temperatures. This may have implications 
for computing in the quantum regime pO| ] . 

For spins coupled to a bath the quantum nature expresses itself in off-diagonal elements of the density matrix. 
These decay after the time T2, which can range up to seconds. In this regime related work extraction setups are 
possible. 

Our results can be phrased in the statement that Maxwell's demon exists: it is the property of quantum entanglement 
in quantum mechanics |]5T|| . They may further have implications for thermodynamics in high-energy physics and the 
early Universe. 

The aim of the paper has been to show that violations of the second law have a natural place in the physics of 
quantum particles that are non- weakly coupled to quantum baths. In this domain we have given conditions for the 
realization of the most notorious objects in the history of physics: perpetuum mobile of the second kind. 



ACKNOWLEDGMENTS 



We thank R. Balian, H. van Beijeren, H. Knops, R. Lipowsky, C. Pombo, H. Spohn, L. Suttorp and M. Wagner 
for discussion and R. Lipowsky also for hospitality at the Max Planck Institute for Colloids and Interfaces in Golm, 
where a part of this work was done. We stress that the responsibility for the presented results and interpretations lies 
with the authors alone. The research of A.E. A. was supported by FOM (The Netherlands) and by NATO. 



APPENDIX A: SMOOTH CHANGES OF THE SPRING CONSTANT 



In this appendix we derive t he wo rk for continuous ch ange s of the spring constant. Hereto we first notice that 
perturbative expression ( 5.55 ), ( 5.57 ) of the exact result ( 5.54 ) can be derived directly by perturbation theory. Let 
us first denote a by ai and expand also to first order in oi ~ ao- We may use 



sm e 



»^o-»0 ^ sin + _ ao)e-'^« ^ = sin q^o - {ai - ao) sin^ </,o e"'*" ^-^T^ 

da 7! V 



Now inserting ( 5.52 ) and taking the large F limit of /(w) from ( 5.44 ), we get from Eq. ( ^.57 ) 



ao — ai / 1 



1 - e 



UJ2 + W 



(A.l) 



(A.2) 



Now let us remember that for t < the spring constant was ao, while for i > it is ai. Thus one could write in 
( ^.36 ): a a{t) = oq + (ai — ao)6'(t). By tre ating the term (ao — ai)9{t)x{t) together with r]{t), one can read off the 
formal solution from the analog of Eq. (5.41), and solve it perturbatively to first order in (ai — ao). It can be verified 
that the result coincides with ( |A.2| ). 

This first order perturbation theory can immediately be generalized for many steps. 



a{t) = afe = (1 - ak)a, {tk < t < tk+i) 



(A.3) 



where to — —00 and ti was taken equal to zero so far, but can be arbitrary. One writes a{t) = ao + X]/c>i('**: 
ak~i)6{t - tk) and gets 



59 



I3{iy) = sin 00 M e' 



1-E 



fc>i 



7W 



0{t-tk) 



I _ g-(wi+ii^)(t-tfc) 2 — g-(cJ2+ii')(t-tfc) 



UJ2 + W 



(A.4) 



At a given instant of time the sum has a finite number of terms because of the 6 functions. 

InEq. ( |A^ ) we considered that effect of many small changes in the spring constant. When we make many changes 
with small a{t) — ak in the domain tk < t < t/c+i, we get 



a 7r7 



2m 



(A.5) 



where ti = ti is the moment of the first change, taken equal to zero so far. Let us write ak = a{Tk) assume and that 
the changes are small. Then the sum can be replaced by an integral, 
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(A.6) 



where the upper integration border could be put equal to oo, since a' vanishes for times less than ti. 
The work needed to make the change Uk-i —> ak at time tk is equal to {uk — ak~^i){x'^{tk)), so it equals 



dWk = (afc-i - ak) { l;T,{a) + ^ 



akCxiO) - ^("j - aj^i)Cxi 



2m ' 



(A.7) 



If there are many steps with small increments, we can go to a continuum limit. Replacing the sum by an integral, we 
obtain the rate of work added to the system 
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dt At 
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where 
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7 da(T) 



di 2m dr 

is the adiabatic (recoverable) rate of work and 

dn ha da(T) 
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is the rate of energy dispersion. 



APPENDIX B: MODERATE CUTOFF FREQUENCY AND FINITE CHANGE OF SPRING CONSTANT 



In this section we address the vanishing of work dispersion at T = 0, without making the approximation of large 
Debye frequency F. Then the full equation ( 5.56 ) has to be employed, rather than the approximation ( 5.57 ). With 



help of an algebraic manipulation program we have checked that at T = the important findings C^T' = 0, see Eq. 



( |6.47 ), and Cp = Cj^ = (c.f. Eq. ( |6.53 )) remain vahd then. So a negative rate of energy dispersion occurs also 



for a finite cutoff. 

Let us mention, however, that the effect is weakened when ao, the amplitude of the change of the spring constant, 
is not very small. This probably affects the maximal number of work extraction cycles. 

If one changes the mass and not the spring constant, the system does not exhibit this interesting behavior, since 



the analog is C, 
low T. 



(0) 



does not vanish then, implying that the leading term in the energy dispersion does not vanish at 
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